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Abstract. Continuing a line of investigation initiated in [11] exploring the 
connections between Jost and Evans functions and (modified) Fredholm de- 
terminants of Birman-Schwinger type integral operators, we here examine the 
stability index, or sign of the first nonvanishing derivative at frequency zero of 
the characteristic determinant, an object that has found considerable use in the 
study by Evans function techniques of stability of standing and traveling wave 
solutions of partial differential equations (PDE) in one dimension. This leads 
us to the derivation of general perturbation expansions for analytically-varying 
modified Fredholm determinants of abstract operators. Our main conclusion, 
similarly in the analysis of the determinant itself, is that the derivative of the 
characteristic Fredholm determinant may be efficiently computed from first 
principles for integral operators with semi-separable integral kernels, which 
include in particular the general one-dimensional case, and for sums thereof, 
which latter possibility appears to offer applications in the multi-dimensional 
case. 

A second main result is to show that the multi-dimensional characteris- 
tic Fredholm determinant is the renormalized limit of a sequence of Evans 
functions defined in [23] on successive Galerkin subspaces, giving a natural 
extension of the one-dimensional results of [11] and answering a question of 
[27] whether this sequence might possibly converge (in general, no, but with 
renormalization, yes). Convergence is useful in practice for numerical error 
control and acceleration. 



1. Introduction 

A problem of general interest is to determine the spectrum of a general variable- 
coefficient linear differential operator L = X)h=o '^a(^)^"' '^ai.x) G R"^", x e W^, 
with prescribed behavior of the coefficients as \x\ — > oo. This arises naturally, for 
example, in the study of traveling- or standing- wave solutions of nonlinear PDEs in 
a wide variety of applications, as described, for example, in the survey articles [35], 
[47], and the references therein. In the one-dimensional case, d = f , a very useful 
and general tool for this purpose is the Evans function [1], [6]-[9], [30], defined as 
a Wronskian £{z) of bases of the set of solutions of the associated eigenvalue 
ODE {L — A)^ = decaying at a; = +oo and x ~ — oo, respectively, whose zeros 
correspond in location and multiplicity with the eigenvalues of L. 
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Among the many applications of the Evans function, perhaps the simplest and 
most general is the computation of the stability index 

r = sgn(4£(0))sgn(£(+oo)), (1.1) 

whose sign determines the parity of the number of unstable eigenvalues A, or eigen- 
values with positive real part Re(A) > 0, where d'^S{0) is the first nonvanishing 
derivative of £{z) at 2; = 0; see, for example, [9], [30], [47]. (A standard property 
of the Evans function is that it may be constructed so as to respect complex con- 
jugation; in particular, it may be taken real-valued for z G M.) A problem that 
has received considerable recent interest^ is to extend the Evans function, and in 
particular the stability index, to the more general setting of multi-dimensions in 
a way that is useful for practical computations. Here, we refer mainly to numeri- 
cal computation, as presumably the only feasible way to treat large-scale problems 
associated with multi-dimensions. 

Various different constructions have been suggested toward this end; see, for 
example, [3], [4], [23]. However, only one of these, the Galerkin approximation 
method of [23] (described in Section 4), seems in principle computable, and the 
computations involved appear quite numerically intensive. (So far, no such com- 
putations have satisfactorily been carried out, though, see the proposed methods 
discussed in [18], [27].) It is therefore highly desirable to explore other directions 
that may be more computationally efficient. 

Here, we follow a very natural direction first proposed in [11]. Specifically, it 
is shown in [11] for a quite general class of one-dimensional operators L that the 
Evans function, appropriately normalized, agrees with a (modified) characteristic 
Fredholm determinant, thus generalizing the classical relation known for the Jost 
function associated with Schrodinger operators of mathematical physics; for further 
discussion of this problem and its history, see [11] and references therein. A central 
point of the analysis is the observation [13] that semi-separability of the integral ker- 
nel of the resolvent operator is the key property of one-dimensional operators that 
makes possible the reduction of the infinite-dimensional characteristic determinant 
to a finite-dimensional determinant expressed by the Jost or Evans function. 

The identification of Evans functions and Fredholm determinants yields a natural 
generalization of the Evans function to multi-dimensions, since the definition of 
(modified) Fredholm determinants extends to higher- dimensional problems. What 
is not immediately clear is whether this extension leads to a practically useful, 
computable formulation of either the Evans function or the stability index. For, 
up to now, the main approach to computation of the characteristic determinant 
was to express it as a Jost function or the usual, Wronskian expression for the 
(one-dimensional) Evans function. 

In the present paper, we extend some of the investigations of [11] in two ways. 
First, we derive in Theorem 2.7 a general perturbation formula for analytically- 
varying modified Fredholm determinants, by which we may express the stability 
index as a product of a finite-dimensional minor (Lyapunov-Schmidt decomposi- 
tion) and a finite-rank perturbation of the original characteristic determinant: that 
is, directly in terms of Fredholm determinants, without reference to a Jost or Evans 
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function formulation. This is carried out in Section 2, with the main result given 
in Theorems 2.3 and 2.7. 

Second, we discuss in Theorem 3.8 an important case when the finite-rank per- 
turbation part, like the original characteristic determinant, may be reduced to 
a finite-dimensional determinant whenever the resolvent of the operator L has a 
semi-separable integral kernel, in particular, in the one-dimensional case. This 
is a consequence of the simple observation that a sum of operators with semi- 
separable integral kernels may be expressed as an operator with a matrix-valued 
semi-separable integral kernel and evaluated in the same way; indeed, the analysis 
of [13] on which this reduction is based is actually presented in the more gen- 
eral, matrix-valued setting. We illustrate this procedure in Section 3 by explicit 
computations for the example of a scalar Schrodinger operator that arises in the 
study of stability of standing-wave patterns of one-dimensional reaction-diffusion 
equations, in the process illuminating various relations between Jost functions and 
characteristic Fredholm determinants. 

To explain our main results in Section 3, we recall that a classical formula by 
Jost and Pais equates the Jost function and the Fredholm determinant of a Birman- 
Schwinger-type operator (cf. Section 3 for its discussion and definitions). Since the 
Jost function is the Wronskian of the Jost solutions this result can be viewed as 
a calculation of the Fredholm determinant via the solutions ^± of the homogenous 
Schrodinger equation that are asymptotic to the exponential plane waves. We 
proved a new formula in this spirit, see (3.99), and computed the derivative of 
the Jost function, that is, the derivative of the Fredholm determinant, via some 
solutions ■^i of a nonhomogenous Schrodinger equation (cf. (3.127)) that, in turn, 
are asymptotic to \I/± (cf. (3.100)). We are not aware of any earlier references 
mentioning these solutions V±- 

In addition, we obtain in Section 3 in passing an elementary proof of an inter- 
esting formula derived by Simon [37] for the Jost solutions ^'±(2:, •) in terms of 
Fredholm determinants. 

Of course, the approach of Section 3 applies equally well to the general one- 
dimensional case, yielding in principle a similarly compact formula for the stability 
index obtained entirely through Fredholm determinant manipulations. However, in 
this paper we do not pursue this any further, leaving this topic for future research. 

More generally, this suggests an approach to the multi-dimensional case by a lim- 
iting procedure based on Galerkin approximation, but carried out within the Fred- 
holm determinant framework. More precisely, we propose in place of the standard 
approach of reducing to a large one-dimensional system by Galerkin approximation 
then defining a standard Evans function, to first relate the Evans function and a 
Fredholm determinant, then evaluate the latter by Galerkin approximation/semi- 
separable reduction. This offers the advantage that successive levels of approxima- 
tion are embedded in a hierarchy of convergent problems useful for error control, 
without the need to prescribe appropriate normalizations by hand. 

Our main result in the multi-dimensional case, generalizing the one-dimensional 
results of [11], is that the sequence of approximate 2-modificd Jost functions T2,j 
generated by Galerkin approximation of the 2-modified Fredholm determinant at 
wave number J agrees, up to appropriate normalization, with the Evans functions 
Ej constructed in [23] for the sequence of onc-dimensional equations obtained by 
Galerkin approximation at the same wave number; see Theorem 4.15. This includes 
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the information that the sequence £j, introduced in [23] as a tool to compute 
a topologically-defined stabiHty index, in fact determines a well-defined 2-modified 
Jost function; that is, an appropriate renormalization ^2, j = e^^Sj of the sequence 
£j converges to a limiting 2-modified Jost function J^2 (cf. Theorem 4.9), a fact 
that is not apparent from the construction of [23]. The Jost and Evans functions 
of course carry considerably more information than the stability index alone. 

We obtain at the same time a slightly different algorithm for computing the 
stability index, which avoids some logistical difficulties of the existing Galerkin 
schemes. We discuss these issues in Section 4, illustrating with respect to the basic 
multi-dimensional examples of flow in an infinite cylinder and solutions with radial 
limits. 

Finally, we note that the ODE systems arising in computation of the Galcrkin- 
based Evans functions £j become extremely stiff as J ^ 00, featuring growth/decay 
modes of order ± J. Thus, numerical conditioning becomes a crucial consideration 
for the large-scale systems that result in multi-dimensions (J ~ 100, as described 
in [18]). It may well be that, for sufficiently large J, direct computation of the 
Fredholm determinant by discretization may be more efficient than either Evans 
function computations or simple discretization of the linearized operator L; see 
Subsection 4.1.7. Efficient numerical realization of this approach would be a very 
interesting direction for future investigation. 

Plan of the paper. In Section 2, we provide a general perturbation formula for 
analytically varying (modified) Fredholm determinants. In Section 3, we use this 
result together with the reduction method of [13] to compute the stability index in 
the case of a one-dimensional real Schrodinger operator. Finally, in Section 4, we 
describe extensions to multi-dimensions. 

2. A GENERAL PERTURBATION EXPANSION FOR FREDHOLM DETERMINANTS 

In this section we describe the analytic behavior of Fredholm determinants 
Aetu{I — and modified Fredholm determinants det2,-H(/ — A{z)) in a neigh- 

borhood of z = with A{ ■ ) analytic in a neighborhood of z = in trace norm, 
respectively, Hilbert-Schmidt norm. Special emphasis will be put on the case where 
[I — A{Q)\ is not boundedly invertible in the Hilbert space H. 

In the first part of this section we suppose that all relevant operators belong to 
the trace class and consider the associated Fredholm determinants. In the second 
part we consider 2-modified Fredholm determinants in the case where the relevant 
operators are Hilbert-Schmidt operators. 

2.1. Trace class operators. In the course of the proof of our first result we re- 
peatedly will have to use some of the standard properties of determinants, such 

as, 



dein{{In-A){In-B)) = deiH{In-A)dein{In~B), A,BeBi{H), (2.1) 
detn'il-H' - AB) = detniln - BA) for all A e B{n, W), B e B{n\ H) (2.2) 

such that BA e Bi(W), AB e Bi{n'), 



and 




(2.3) 
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Since 

'«-^=(o j:X)'{'o iJd^C^ t)- ("> 

Finally, assuming A,Bg Bi{H), we also mention the following estimates: 

|det„(7„ - A)\ < exp(|m|B,(„)), (2.5) 
|det„(7„ -A)- detniln - B)\ < \\A - (2.6) 

xexp(||A||B^(„) + ||B||B^(„)+l). 

Here H and H' arc complex separable Hilbert spaces, B{Tl) denotes the set of 
bounded linear operators on H, BpiTi), > 1, denote the usual trace ideals of 
B{H), and I-h denotes the identity operator in H (similarly, Ik abbreviates the 
identity operator in C'^). The ideal of compact operators on Ti. will be denoted by 
Boa{H). Moreover, detp^7^(/7^ — A), A £ Bp{H), denotes the (p-modified) Fredholm 
determinant of I-^ — A with deti^-^(/-^ — A) = det-^(/-^ ~ ^)) ^ € Bi{H), the 
standard Fredholm determinant of a trace class operator, and tT-niA), A G Z?i(7i), 
the trace of a trace class operator in H. Finally, + in (2.3) denotes a direct but 
not necessary orthogonal direct decomposition of H into K. and the fc-dimensional 
subspace C'^. These results can be found, for instance, in [14], [17, Sect. IV. 1], [32, 
Ch. 17], [36], [38, Ch. 3]. In the following, <t{T) denotes the spectrum of a densely 
defined, closed linear operator T in H, and <Jd{T) denotes the discrete spectrum of 
T (i.e., isolated eigenvalues of T of finite algebraic multiplicity). 

For the general theory of (modified) Fredholm determinants we refer, for in- 
stance, to [5, Sect. XI.9], [14], [15], [16, Ch. X.III], [17, Ch. IV], [36], and [38, Sects. 
3, 9]. 

Hypothesis 2.1. Suppose A{ ■) E BiCH) is a family of trace class operators on Ti 
analytic on an open neighborhood flo G C of z = in trace class norm \\ ■ WbiCH)- 

Given Hypothesis 2.1 we write 



A{z) = Ao + Aiz + O(z^) for zG^o suSiciently smah, Ag e BUH), £ = 0,1 

(2.7) 



We start by noting the following well-known result. 
Lemma 2.2. Assume Hypothesis 2.1 and suppose [lu — ^o)~^ G ^^)- Then, 

detn {Ih - A{z)) =^ det„ {In -Aq)- detn {In - Aq) tr« ((/„ - Aq) Ai) ^ + O {z^) . 

(2.8) 

Proof. This follows from 
det„(7„ - A{z)) detn{{I-H - Ao)[ln - {In - Ao)-^A,z + 0{z^)]) 

=^ detn{In - ^o) [l - tr«((/„ - Ao)-^A^)z + 0{z^)] , ^^'^^ 
where we used the fact that 



detn{In - Bz) = exp 



V tr^(^') ,fc 



for |z| sufficiently small (2.10) 



km 

withBGBi(W). □ 
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Next we turn to the case where {In — Aq) is not boundedly invertible. Before 
we state the analog of Lemma 2.2 in this more general setting, we need some 
preparations. 

We temporarily assume 

Ao e B{n) and 1 G ad{Ao), (2.11) 

and abbreviate by Pq the Riesz projection associated with Aq and the discrete 
eigenvalue 1 of Aq, 

Po = ^^£ dC{Ao-CIn)-\ (2.12) 

where Co denotes a sufficiently small counterclockwise oriented circle centered at 1 
such that no part of f7(^o)\{l} intersects Co and its open interior. We denote by 

no = dim(ran(Po)) (2.13) 

the algebraic multiplicity of the eigenvalue 1 of ^o- In addition, we introduce the 
quasinilpotent operator Dq associated with Aq and its discrete eigenvalue 1 by 

Do = iAo-In)Po (2.14) 

such that 

Do = DoPo = PoDa = PoDqPo- (2.15) 
In the following we denote by A(zo;ro) C C the open disc centered a,t zq & C 
of radius ro > and hy C{zo;ro) = dA{zo;ro) the counterclockwise oriented circle 
of radius tq > centered at zq. Assuming that the analytic family of operators 
A{-) satisfies Hypothesis 2.1, we prescribe an eq > and choose a sufficiently small 
open neighborhood Oq of ^ = such that all eigenvalues Xj{z) of A{z) for z €ilo, 
which satisfy Aj(0) = 1, 1 < j < fo for some I'o G N with I'o < no, stay in the disc 
A(l; eo/2). Moreover we assume that Qq is chosen sufficiently small that no other 
eigenvalue branches of A(z), z G Qq, intersect the larger disc A(l;eo). Introducing 
the Riesz projection P{z) associated with A{z), z G Qq (cf., e.g, [21, Sect. III. 6]), 

P{^) = ^i dC{A{z)-(:in)-\ zeQo, (2.16) 

27ri /c(l;eo) 

then P( • ) is analytic in flo and we expand 

P{z) =^ Po + P^z + 0{z^) for 1^1 sufficiently small. (2.17) 

Moreover, we introduce the projections 

Q{z) = In - P{z), zeno, Qo = Iu-Po, (2.18) 

and expand 

Q(z) =^ Qa + Qiz + 0{z'^) for \z\ sufficiently small. (2.19) 

Since Piz)"^ = P{z), (2.17) implies 

PoPi + PiPo = Pi and hence PoPiPo = 0. (2.20) 
Following Wolf [45] wc now introduce the transformation 

Tiz) = PoP{z) + QoQ{z) = PoP{z) + [In-Po][In-P{z)], z^Qo, (2.21) 
such that 

PoT(^) = T(^)P(^), QoT{z) = T{z)Q{z), zgQq. (2.22) 
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In addition, for \z\ sufficiently small, 

T{z) = In + iPoPi-PiPo)z + 0{z^), (2.23) 

T{z)-^ ^In-iPoPi-PiPo)z + 0{z^), (2.24) 

and hence 

Po = T{z)P{z)T{z)-\ Qo = T{z)Q{z)T{z)-' (2.25) 
for \z\ sufficiently small. Below wc will use (2.25) to reduce determinants in the 
Hilbert space P{z)'H to that in the fixed Hilbert space PoH. 

Next, we introduce the following notation: We denote by 5 e S(PoW, C"°) the 
boundedly invertible linear operator which puts the nilpotent operator PqDqPq 
into its no x no Jordan canonical form SPoDoPoS~^ , and abbreviate by i^o the 
number of entries 1 in the canonical Jordan representation of SPoDoPoS~^, where^ 
< 1^0 < nQ — I. Moreover, wc denote by Ai the (no — vq) x (no — z/o)-matrix 
obtained from S'Po^i-Po'S'"^ by striking from it the i^o columns and rows in which 
SPoDoPoS~^ contains an entry 1. With this notation in mind, we now formulate 
our first abstract result on expansions of Fredholm determinants detniln — ^(-2)) 
as z ^ 0, with [I-h — A{0)] not boundedly invertible in H: 

Theorem 2.3. Assume Hypothesis 2.1 and let 1 G ad{Ao). Then, given the nota- 
tion in the paragraph preceding this theorem, 

detniln - A{z)) = [detQ^nihoH - QoAoQo) + 0{z)] 

X (-l)"''detp„„(Po£'oJ'o + Po[Ai + 0{z)]Poz) 

detgoW^QoW - QoAoQo)detcno-.o (Ii)(-^)"°-"'' + 0(0"°-^°+!). (2.26) 

In the special case where 1 is a semisimple eigenvalue of Aq {i.e., where Dq = 
and uq = 0) one obtains, 

det„(/w - A{z)) =^ detQMlQon - QoAoQo)detpMPoAiPo){-zr° + 0(z"«+i) 
=^ detniln - Po - ^o)detp„„(PoAiPo)^"'' + 0{z^°+^). 

z— >0 ^ ' 

(2.27) 

In particular, ifl is a simple eigenvalue of Aq {i.e., if no = 1, Dq = 0, and vq = 0), 
one obtains 

detn{lH - A{z)) detn{lH - f - ^o) detp„w(f'oAifb)-2 + 0{z'^). (2.28) 
Proof. Since 

H = P{z)H + Q{z)n, P{z)A{z) = A{z)P{z), P{z)Q{z) = Q{z)P{z) = 0, 

(2.29) 

one computes using (2.22), 

det„(JH - A{z)) = det„(/„ - T{z)A{z)T{z)-^) 

= detn{lH - T{z)[P{z)A{z)P{z) + Q{z)A{z)Q{z)]T{z)-') 
= detn{l-H - PoT{z)A{z)T{z)-^Po - QoT{z)A{z)T{z)-^Qo) 

^In particular, ua equals the sum of the dimensions of all nontrivial (i.e., nondiagonal) Jordan 
blocks in the canonical Jordan representation SPoDoPqS^^ of PoDoPq. Thus, SPoDoPoS~^ 
contains i^o entries 1 at certain places right above the main diagonal and O's everywhere else 
including on the main diagonal. 
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= detp,H{lp„n - PoT{z)A{z)T{z)-^Po) 

X detQ„H{lQoH - QoT{z)A{z)T{z)-^Qo). (2.30) 
Using (2.20), (2.7), (2.23), and (2.24), one computes 

PoT{z)A{z)T{z)-^Po = PoAoPo + PoAiPoz + PoO(z^)Po, (2.31) 

QoT{z)A{z)T{z)-'Qo = Qo^oQo + QoO(^)<3o (2.32) 

for \z\ sufficiently small. Relation (2.32) implies 

detgMlQon - QqT{z)A{z)T{z)-^Qo) = detQ^n{lQ^H - Oo^oOo) + 0{z) ^ 

«— »o 

(2.33) 

for \z\ sufficiently small, and hence we next focus on the first factor on the right- 
hand side of (2.30). Applying (2.14), (2.15), and (2.31) one obtains 

Aeip,H{lpon - PoTiz)A{z)T{z)-'Po) 

=^ dctpMlPoH - PoAoPo - Po[Ai + 0{z)]Poz) 

= i-ir^detpMPoDoPo + Po[Ai + 0{z)]Poz). (2.34) 

Next, let 5 e Bi^PoHjC"") be the transformation which puts PqDqPo into its 
Jordan canonical form Dq = SPoDoPoS~^ and denote 

Ai{z) = SPq[Ai+0{z)]PoS-^ = Ai+0{z). (2.35) 

z— ►O z— >0 

Then, 

detp„„(7p,„ - PoT{z)A{z)T{z)-^Po) 

=^ {-ir°detpMPoDoPo + Po[Ai + 0{z)]Poz) 

z— »0 

= (-l)"Metc"o(5o + Ai(0)^) 

(-l)"Metc"o {Do + Aiz + 0{z')) 
= (-z)"»-^°detc»o-o(^i + 0(z)) 

= (-2)"°-''«detc.o-o(^i) +0(^"°"""+'), (2.36) 

z— >0 \ / \ 

by applying the Laplace determinant expansion formula (cf., e.g., [43, Sect. 3.3]) to 
detc"o {Dq + A\z + 0(2;^)) with respect to the columns in £)q which contain a 
1. Combining (2.36) and (2.33) then proves (2.26). 

If 1 is a scmisimplc eigenvalue of A^ and hence Dq = 0, ;/o = 0, the first line on 
the right-hand side in (2.27) is clear from (2.26). To prove the second line in the 
right-hand side of (2.27), we recall that 

Po^o^b = -Po^o = -4oPo = -Po and Qo^oQo = Qo^o = {In - Po)Ao = Aq ~ Pq, 

(2.37) 

and hence, 

detn{lH -Po- Ao) = detniln - Po - PoAqPo - QoAqQo) 
= det-H{-Po + Qo- Qo^oQo) 

= i-ir°detQ,H{lQon - QoAoQo) + 0. (2.38) 
The special case no = 1 in (2.27) then yields (2.28). □ 
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Remark 2.4. Since i^o can take on any particular value from to no — 1 and Ai 
is generally independent of Aq and hence Do, the power no — fo of z in (2.26) can 

take on any value between 1 and uq. 

2.2. Hilbert— Schmidt operators. Next, wc treat the case of 2-modificd Fred- 
holm determinants, where all relevant operators are only assumed to lie in the 
Hilbert-Schmidt class. In addition to (2.1)-(2.3) we recall the following standard 
facts for 2-modificd Fredholm determinants dct2.n{I'H ~ A), A € B2{'H) (cf., e.g., 
[15], [16, Ch. XIII], [17, Sect. IV.2], [36], [38, Ch.' 3]), 

dct2,n{lH -A)= detniilH - A) exp(A)), A e S2(H), (2.39) 

det2,w((Jw - A){In - B)) = det^M^n - A) det2,w(Jw - B) e-t'^«(^^), (2.40) 

A,BeB2{n), 

det2,H{In -A) = detniln - A) e*'-"^^), A e Si(H), (2.41) 
det2,n'{In' - AB) = det2,n{In - BA) for all A e B{n,n'), B e B{H',n) 

such that BA e B2{n), AB e S2(W)- (2-42) 
Moreover, in analogy to (2.10) one now has 



det2,H{lH - Bz) = exp 



fe=2 ^ 



-Z^ 



for \z\ sufficiently small, B e B2{'H). 



(2.43) 

Finally, assuming A,Bg B2{'H), we mention some estimates to be useful in Section 
4: For some C > 0, 

|det2,^,(/w - A)\ < exp {C\\A\\l^^^^) , (2.44) 
|det2,H(Jw -A)- det2,H(/w -B)\<\\A- By^in) (2.45) 

xexp(C[||A||B,(„) + ||B||B,(„) + l]2) 

Hypothesis 2.5. Suppose A{-) € -62 (H) is a family of Hilbert-Schmidt operators 
on H analytic on an open neighborhood fio C C 0/ .2 = in the Hilbert-Schmidt 
norm || ■ We^iH)- 

Given Hypothesis 2.5, we write again 

A(z) = ^0 + Aiz + 0(z'^) for zgQo sufficiently small, Ae e B2(W), ^ = 0, 1. 

z— »0 

(2.46) 

Wc start with the analog of Lemma 2.2. 

Lemma 2.6. Assume Hypothesis 2.5 and suppose (/-^ — ^0)^^ G B{Ti.). Then, 

det2,H(lH - A(z)) = dct2,w(/w - Aq) 

-det2,Hi.In-Ao)tin{i.In-Ao)-^AoAi)z + 0{z^). 

(2.47) 

Proof. This follows most easily from rewriting (2.8) in terms of modified Fredholm 

determinants det2,H{lH+ • ), using (2.41), and then approximating Hilbert-Schmidt 
operators by trace class (or finite-rank) operators (cf., e.g., [17, Theorem HI. 7.1]). 
Indeed, one computes using (2.41) repeatedly, 

det2,w(/w - A{z)) = detniln - A(^))e*--«(^(^)) 
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= detn{{In - Ao)[ln - (In - Ao)-\A{z) - Ao)])e''-«(^(^)) 
= det«(/« - Ao)det„(/„ - {In - A^)-\A[z) - Ao))e*'-«(^(^» 
= det2,H(/w - Ao)e-*'"«(^°)det2,w(/H - (/« - A^Y^^A^z) - A^)) 

= det2,H(-fw - Ao)det2,H(/w - ijn - Aq)-^{A{z) - A^)) 

= det2,w(/w - Ao)det2,H{ln - {In - Ao)-\A{z) - Aq)) 

= det2,w(/H - Ao)det2,w(/w - (/« - Aoy^Aiz + 0(z'^)) 

X g-tr„((7„-Ao)-iAoAiz+0(z=)) 

det2,w(/-H - Ao)dct2,w(/-H - (/« - Ao)"^AoAiz + 0{z^)) 

X [1 - tr„((/„ - Ao)-^AoAi)z + 0{z^)] 
= det2, n{In - Aq) 

z— »0 

- det2,^^(/w - Ao)tr„((/„ - Aor'AoAi)z + 0{z^). (2.48) 
Here we used (cf. (2.43)) 

det2,w(/w - B{z)z) 1 + 0{z^) (2.49) 

for B{ ■ ) analjftic in S2(W)-norni near z = 0. □ 

In exactly the same manner one obtains the Hilbert-Schmidt operator version of 
Theorem 2.3. Again we rely on the notation introduced in the paragraph preceding 
Theorem 2.3. 

Theorem 2.7. Assume Hypothesis 2.5 and let 1 G adiAo). Then, 
det2,w(/w - A{z)) [dct2,QMlQon - QaAoQo) + 0(;j)]e"« 

X i-ir'dctpMPoDQPo + Po[Ai + 0{z)]Poz) (2.50) 
det2,Qow(/QoH - Qo^oQo)e"Metcno-^o (Ii) (-z)"""''" + 0(2"°--°+i) . 

/n special case where 1 is a semisimple eigenvalue of Aq {i.e., where Dq = 

and pq = 0) one obtains, 



det2,H(/-H - A{z)) = det2^Q„n{lQoH - QoAoQo)e''°detpMPoAiPo){-z) 



no 



z->0 



+ 0(z"°+i) (2.51) 
= det2,H(/H -Po- Ao)e"°detp„„(Po^iPo)^"° + 0(^"''+'). 



z->0 



In particular, ifl is a simple eigenvalue ofAo {i.e., if no = 1, Dq = 0, and uq = 0), 

one obtains 

det2,H{In - A{z)) = det2,H{In - Po - Ao)detp,n{PoAiPo)z + 0{z^). (2.52) 

z— >0 
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Proof. Again this follows from rewriting (2.26)-(2.28) in terms of modified Fred- 
holm determinants det2,n{In + using (2.41), 

tr„(^(^)) = trH((Po + Qo)T{z)A{z)T{z)-\Po + Qo)) 

= tTH{PoT{z)A{z)T{z)-^Po) +tvn{QoT{z)A{z)T{z)-^Qo) 

= tvp^n{PoTiz)Aiz)Tiz)-'Po) + tvQ^n{QoTiz)Aiz)Tiz)-'Qo) 

= tvp^-n{PoAoPo) + 0{z)+tvQ,n{QoT{z)A{z)T{z)-^Qo) 

=^ no + 0{z) + tvQ^n{QoT{z)A{z)T{z)-'Qo) , (2.53) 

and then approximating Hilbert-Schmidt operators by trace class (or finite-rank) 
operators (cf., e.g., [17, Theorem III. 7.1]). Here the fact that trp^niPoAoPo) = 
no follows from the canonical Jordan structure of SPoAoPoS~^ . Explicitly, one 
computes 

det2,w(/w - A{z)) = detniln - A(^))e*--«(^(^)) 
= detp„„(/p„„ - PoT{z)A{z)T{z)-^Po) 
X det2,Q,H{lQon - QoT{z)A{z)T{z)-^Qo) 

X e~^'"='o'n{QoT(z)A(z)T(z)-^Q„)^trT^{A{z)) 

= detp,n{lpon - PoT{z)A{z)T{z)-'Po)[det2,QMlQon - QoAoQo) + 0{z)] 

y ^no+0(z) 

= (-l)"°detp„„(Po£'oPo + Po[Ai + 0{z)]Poz)) 

z— >0 

X [det2,Qow(/QoW - QoAoQo) + 0(^)]e"°. (2.54) 

□ 

2.3. Higher modified determinants. It is now obvious how to proceed in con- 
nection with higher modified Fredholm determinants detp( • ), p G N, n > 3, when 
the family A{-) in Hypothesis 2.5 is replaced by one analytic near 2; = in Bp{H)- 
norm. We omit further details at this point. 

3. One-Dimensional Reaction-Diffusion Equations and Schrodinger 

Operators with Nontrivial Spatial Asymptotics 

As an elementary illustration of formula (2.28) and a view toward applications 
to reaction-diffusion equations, we now illustrate the abstract results of Section 2 in 
the context of one-dimensional Schrodinger operators. (For background literature 
on reaction-diffusion equations we refer, e.g., to [29], [33], [39].) 

To motivate our considerations of one-dimensional Schrodinger operators in this 
section, we start with a brief discussion of a simple model for one-dimensional scalar 
reaction-diffusion equations of the type 

m = Wxx + f{w), t>0, w{0) = Wq. (3.1) 

Here, we assume, for simplicity, 

/ e Ci(R), /,/' G L°°(R), Wo e i^^(M) nC°°(]R), e ff2(M), (3.2) 

with w = w{x,t), {x,t) e R X [0,T], for some T = T{wq) > 0, a mild solution of 
(3.1) satisfying w e Cb{[0,T];L'^{R;dx)). 
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As usual, Cb{[0,T];L^{R;dx)) denotes the space of bounded continuous maps 
from [0, T] with values in L^(IR; dx) and the sup norm 

\\v\\c,i[o,T];L^R;dx)) = sup \\v{s)\\L2(R-dx), V G Cbi[0,T]; L"^ {R; dx)). (3.3) 

se[o,T] 

We refer, for instance, to [24] for more details in this context (see also, [19, Sect. 
3.2]). Also, H'^{R), m e N, denotes the standard Sobolev spaces of regular distri- 
butions which together with their derivatives up to the mth-order lie in _L^(R; dx). 

Assuming that U is a, stationary (steady state) solution of (3.1), that is, Ut = 0, 
and hence 

U" + f{U) = 0, (3.4) 
we now linearize (3.1) around U and obtain the linearized problem, 

vt = Lv, (3.5) 

where L in L^(M; dx) is given by 

L=-^ + f'{U), dom{L) = H\R), (3.6) 

Since U' G H^{R) n C°°(IR), one infers upon differentiating (3.4) with respect to x 
that 

L{U') = iU')" + f{U)U' = [U" + /([/)]' = 0, (3.7) 

and hence, G crp(L). Here crp( • ) abbreviates the point spectrum (i.e., the set of 
eigenvalues). Thus, the Schrodinger operator 

H = -L = -^-f'{U), dom{H)=H\R), (3.8) 

in (M; dx) with the potential 

V{x) = -f\U{x)), xGR, (3.9) 

has the special property of a zero eigenvalue, that is, 

e ap{H), (3.10) 

and we will be studying the situation where in addition is a discrete (and hence 
simple) eigenvalue of H in great detail in the remainder of this section. In particu- 
lar, under appropriate assumptions on the the "potential" term —f'{U) in H, the 
Birman-Schwinger-type integral operator K{z) associated with H and the complex 
eigenvalue parameter z £ C will be a trace class operator with K{0) having the 
eigenvalue 1. Consequently, the behavior of detj;,2(R.(ij,)(Jj;,2(R.(ij,) — K{z)) for z in 
a sufficiently small open neighborhood of z = 0, which determines linear stability 
of (3.1) around its stationary solution U (cf. [19, Sect. 5.1, Ex. 6, Sect. 6.2], [25, 
Sects. 9.1.4, 9.1.5]), can be determined in accordance with Theorem 2.3 or alterna- 
tively, directly from Jost (respectively, Evans) function considerations, as discussed 
in Remark 3.5. 

We conclude these illustrations with the following elementary example: 
Example 3.1. Let n G N, n > 2, c > 0, k > 0, a; G M, and consider 

f„{w) = -{n - l)^K^w + {n- l)nK2c-2/("-i)u,("+i)/("-i), (3.11) 
Unix) = c[cosh(Kx)]-"+^ (3.12) 
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Then the potential Vn{x) = — f {Un{x)) in the corresponding Schrddinger operator 
Hn = —{d'^/dx^) + Vn coincides with a particular family of n-soliton Korteweg-de 
Vries (KdV) potentials 

Vn{x) = -f'{Un{x)) = (n- l)^K^ - n{n + l)K2[cosh(Kx)]-^ (3.13) 

lim Vn{x) = {n- I^k^ > 0, (3.14) 

X— >±oo 

and the zero-energy eigenfunction U'^ satisfying Hn{U^) = is given by 

U^{x) — —{n — l)Kcsinh(Kx)[cosh(Kx)]~". (3.15) 

It seems a curious coincidence that Vn should coincide with a particular family 
of n-soliton KdV potentials (there are many other such n-soliton KdV potentials, 
cf., e.g., [10, Example 1.31]) in this reaction-diffusion equation context. 

Motivated by these considerations, we now start to investigate one-dimensional 
Schrodinger operators with a scalar potential displaying a nonzero asymptotics as 
\x\ cx). 

Hypothesis 3.2. Assume V: R^M. is Lebesgue measurable and satisfies 

dx \V{x) — Vool < oo for some V^o > 0. (3.16) 



Assuming Hypothesis 3.2 we introduce 

W{x) = V{x) - Foo, (3.17) 

W{x) = vix)u{x), u{x) = sgn{W{x))v{x), v{x) = \W{x)\^/'^, (3.18) 

for a.e. x € R, and define H to be the (maximally defined) self-adjoint realization 
in L^(R; dx) of the differential expression r = —{d^/dx^) + y(x), x £ K, obtained 
by the method of quadratic forms, or cquivalcntly, by using the limit point theory 
for sclf-adjoint 2nd order ordinary differential operators, 

Hf = Tf, (3.19) 

/ e Aom{H) = {g G L^{R;dx) \ g,g' G ACioc(M);Tg £ L^{R;dx)}. (3.20) 

(Here / denotes differentiation with respect to x G R.) This also implies 

dom(iJi/2) = i7i(M). (3.21) 

We also introduce the self-adjoint operator in L'^{R;dx) associated with the 
differential expression r'^^ = —{d'^/dx^) + Voo, x G K, replacing y( ■ ) by its asymp- 
totic value Voo (in the sense of (3.16)), 

^(o)y^^(o)j^ / e dom(i?(°)) = i?2(R). (3.22) 
By well-known results, (3.16) implies 

aess(if) = ae.s(i/'°') = [l^oo,oo), adiH) C i-^,V^), (3.23) 

where (Tess( • ) and cr(i( • ) denote the essential and discrete spectrum, respectively. 
Applying the Birman-Schwinger principle (cf., e.g., [12, Sect. 3]), 

H^{Xj) = Aj*(Aj), Xj < Voo, Xj G (Td{H), *(Aj) G dom(if), (3.24) 

is equivalent to 

K{Xj)^Xj) = $(A,), Xj < Voo, ^Xj) G L^{R; dx), (3.25) 
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with equal finite geometric multiplicity of either eigenvalue problem (3.24) and 
(3.25). In particular, in this special one-dimensional context, the eigenvalue Xj of 
H as well as the eigenvalue 1 of K{Xj) are necessarily simple. Here we abbreviated 



K{z) = -u{mo) -zlLHR-d.)) V ^gC\[Ko,oo), (3.26) 

with T denoting the operator closure of T. Wc recall that the integral kernel 
— 2:/l2(r.£;3,)) (x, x') of the resolvent of H^^^ is explicitly given by 

{H^'^->-zI,^^^,,,r'M = (V2)(^-Ko)-^/^e*(^-^-)^^^l--'l, ^^^^^^ 
z gC\[Voo,oo), Im{{z-Voo)^/'^) >0, x,x' gR, 

and hence 

K{z)eBi{L^{R;dx)), z€C\[V^,oo), (3.28) 
since K{z) is the product of the two Hilbert-Schmidt operators (cf. [38, Ch. 4]) 



u{H^"^ - zlL2^t,,d,)) ^'^ and - z/^^dj^.^,)) ^'""v (3.29) 

(no operator closure necessary in the first factor). In addition, 4'(Aj) and ^•(Aj) 
are related by 

$(Aj,a;) = ?i(.T)*(Aj,a;), a; e R, (3.30) 
and we note that ^'(Aj, • ) is also bounded, 

*(Aj, •) G i°°(]R) (3.31) 

(in fact, even exponentially decaying with respect to x by standard iterations of the 
Volterra integral equations (3.43)). Here the discrete eigenvalues (Td{H) = {Xj}j^j 
oiH, with J C N an appropriate (finite or infinite) index set, are ordered as follows: 

Ai < A2 < ••• < 14o. (3.32) 

Moreover, one obtains 



K{zy = -v{H(.^) - zlL^(^,a.)) 'u, ^GC\[yoo,oo), (3.33) 

and hence 

K{zy = SK{z)S-\ z e C\[Voo, oc), (3.34) 
where S denotes the unitary operator of multiplication by sgn(VK( • )) in _L^(]R; dx), 

{Sf){x) = sgn{W{x))f{x) for a.e. x eR, / G L^{R; dx). (3.35) 

Thus, 

KiXjY^iXj) = $(Aj), Xj < Voo, ^Xj) G L'^{R;dx), (3.36) 

where 

^Xj,x) = C{Xj){S^Xj)){x) = C{Xj)v{x)^{Xj,x) for some C{Xj) G C\{0}. 

(3.37) 

To fix the normalization constant C{Xj) wc require 

($(A,),<E>(A,))^,(^^,^^ = 1. (3.38) 

This yields 



C{Xj) 



ljxW{x)\nXj,x)f =(t^*(A,),«*(A,))-.\j,^,,) (3.39) 
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and hence 



/ 

Jr 

j 



dxW{x)\^{Xj,x)[' 



dxW{x)\^{\j,x)[' 



sgn{W{x))^{Xj,x) (3.40) 



v{x)^{\j,x). (3.41) 



In particular, the corresponding one-dimensional Riesz projection P{Xj) onto the 
eigenspace corresponding to the eigenvalue 1 of K{\j) is then of the simple form 



P(A,) = ($(A,),.) 



(3.42) 



3.1. The Jost— Pais derivative formula. Next, consider the generalized Jost- 
type (distributional) solutions 



^±{Z,X) = 6=" 



/■iboo 

- / dx 

J X 



,sm{{z-V^Y'^{x-x')) 



[V{x')-Voo]^±{z,x'), (3.43) 



z e C\[V^, oo), lm((z - V^f'^) > 0, a; G : 

oiL'^{z) = z*(2). Then, 

^±{X,x) are real- valued for A < Vqo, a; G M, 



and 



^± (A, a;) > for A < Voo and ±x sufficiently large. 
The Jost function T associated with H is then given by 

Wr(*_(z),*+(0)) 



Hz) 



2i{z - v^y/^ 



= 1 + 



2(^ - 14o)V2 A 



/ 



dx e 



[V{x)-V^]^±{z,x), 



(3.44) 
(3.45) 

(3.46) 

(3.47) 



2eC\[Ko,w), 

where Wr(/,s()(a;) = f{x)g'{x) - f{x)g{x), x € R, f,g € Ci(IR), denotes the 
Wronskian of / and g. The limits to the real axis 

limJP(A±ie) =T(X±iO) 
exist and are continuous for all A G K\{K>o}- In addition, one verifies 



(3.48) 



I >00 



(3.49) 



Moreover, one can prove the following result originally due to Jost and Pais [20] 
(for Voo — 0) in the context of half-line Schrodinger operators. The actual case 
at hand of Schrodinger operators on the whole real line (again for Voo = 0) was 
discussed by Newton [26] and we refer to [13] for more background and details, 

T{z) = dcti2(R.rf^)(/i2(R.rf^) - K(z)), z G C\[Ko, oo). (3.50) 

Since we are interested especially in the z-derivative of T{z) at a discrete eigen- 
value of we now prove the following result. Lemma 3.3. For the remainder of 
this paper we abbreviate differentiation with respect to the complex-valued spec- 
tral parameter z G C by • (to distinguish it from differentiation with respect to the 
space variable a; G K). 
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Lemma 3.3. Assume Hypothesis 3.2 and z € C\[yoo,c»). Moreover, let Xj < Voo, 
Xj e ad{H). Then, T{Xj) = and 

= oiv ~\ M/2 / dx^_{X„x)^+{X,,x). (3.51) 

Proof. Consider 

■^l{z,x) = [V{x)-z]^±{z,x), -^'^{z^x) = [V{x)-z]^'±{z,x)-^±{z,x) (3.52) 
(we recall that • abbreviates d/dz) to derive the identities 

^m{^-{z,x),^\{z,x)) = -^-{z,x)^+{z,x), (3.53) 

^Wr(**.(z,a;),vE'+(z,a;)) = ^-(z, a;)1'+(z, x). (3.54) 
Then one obtains for all i? > 0, 
Wr(*_(z), ^'+{z)){x) + Wr(*l(0), ^+{z)){x) 

- Wr(*_(z), ^'+{z)){R) - Wr(**_(z), ^+{z)){-R) 

= -J dx' -^Wv(<i!-{z,x'),^\{z,x')) + j dx'^Wr(**_(z,a;'),*+(^,a;')) 

= / rfa;'*_(^,a;')^'+(2,a;'). (3.55) 

Next, wc note that 

£wr(*_(z),*+(z)) = Wr(*_(z),*;(z))+Wr(*'_(z),*+(z)), (3.56) 

and choosing z = Xj < 14o, Xj e a^iH), one concludes ^'±(Aj, ■ ) G L^(M; rfa;) and 
hence, 

^Wr(*_(z),vl/+(z))|^^^, = jrda;M/_(A,-,a:)M/+(A,-,a:). (3.57) 
Together with (3.46) and 

2i(A,- - V^y/^J^{Xj) = Wr(*_(A,),*+(A,)) = 0, (3.58) 
this establishes (3.51). □ 
Next, we specialize to the case Xj = and hence assume 

e ad{H). (3.59) 

In this context we then denote 

*o = *(0), $0 = $(0), $0 = ^(0), Co = C(0), Ko = K{0), Pq = P(0), etc., 

(3.60) 

and without loss of generality (cf. (3.44)) we assume that ^'o(a;) is real-valued for 
all X e R. 

We summarize the results for !F*{0): 

Lemma 3.4. Assume Hypothesis 3.2 and suppose G ad{H). Then, J^{0) = and 
^*(0) = — 172 / dx^-{0,x)<b+{0,x) (3.61) 
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_ I < if \E'±(0, • ) has an even number of zeros on R, 
1 > if *±(0, • ) has an odd number of zeros on M, 

{< if is an odd eigenvalue of H, 
> if is an even eigenvalue of H. 

Here the eigenvalues {Aj jjgj, J C N an appropriate {finite or infinite) index set, 
are ordered in magnitude according to Xi < X2 < ■ ■ ■ < Voo {cf. (3.32)). 

Proof. Equation (3.62) is immediate from (3.45) and ^'+(0, x) = cvl'_(0, x) for some 
c e K\{0} (cf. (3.97)). Relation (3.63) is a direct consequence of the fact that H 
is bounded from below, H > \ilL^(^-dx)^ the discrete eigenvalues of H are in a 
one-to-one correspondence with the zeros of on [Ai.V^oo) (the zeros necessarily 
being all simple), and the fact that J^{z) — > 1, 2; G C\[V'oo,oo) (cf. (3.49)). □ 

Remark 3.5. Since r{\ ± iO) — > 1 by (3.49), = and JP*(0) > imply 

A— *±oo 

that H has at least one negative eigenvalue and hence L = —H (cf. (3.8)) has 
at least one positive eigenvalue. This implies linear instability of the stationary 
solution U in the context of the reaction-diffusion equation (3.1) identifying V{x) 
and —f'{U{x)), x £ R (cf. (3.9)). For, it is easily seen by consideration of the 
standing-wave equation (3.4), a scalar nonlinear oscillator, that the derivative U' 
of a pulse- type solution has precisely one zero, whereas 'ii± are nonzero multiples 
of the zero eigenfunction U'. For discussion of spectral stability and some of its 
applications, we refer, for instance, to [30], [31], [35], [47], [48], and the references 
cited therein. An equivalent formula for J-'* yielding the same conclusions may 
be derived in straightforward fashion by Evans function techniques, following the 
standard approach introduced in [6], [7], [8], [9]. We recall (cf. [11]) that Jost and 
Evans functions, suitably normalized, agree. 

3.2. Fredholm determinant version. We now turn to the connection with the 
abstract approach to the asymptotic behavior of Fredholm determinants presented 
in Section 2. 

Applying (2.28) to (3.50) one then obtains 

^•(0) = deti2(R.(j3,)(7i2(R.rf3,) - Kq- Po) detpgi2(R.rf3,)(Poii'iPo), (3.64) 

where 

K(z) = Ko + Kiz + O(z^) (3.65) 
^— »o 

with 



Ka = -uH(o) \, (3.66) 
Ki = K*{z)\^^^^ = -uH(")~'^v. (3.67) 

3.2.1. Evaluation of the second factor. We start by determining the second factor 
'i^ipoL'^(TS.;dx){PoKiPo) on the right-hand side of (3.64): 

Theorem 3.6. Assume Hypothesis 3.2 and suppose G ad{H). Then, 

detp^L2(^^.da:){PoKlPo) = -W^oWh^^.^x)/ {■"^0,U^o)mR;dx) (3-68) 

= [ll*0lli2(R;da;) +14011*01112(8.^2;)] ^||*0||i2(R;d3;) > 0. 

(3.69) 
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Proof. First, we choose a compactly supported sequence Wn € C^(K), n £ N, such 
that 

W„ = V„-Voo= u„v„, u„ = sgn{W„)vn, v„ = |W„|^/^ n e N, (3.70) 

and 

Hm \\vn - v\\L2(R-^dx) = 0. (3.71) 

Given Wn, wc introduce the self-adjoint operator sum i7„ = H'-^^ +Wn in dx) 

defined on the domain ff^(R) associated with the differential expression L„ = 
— {(P /dx^) + Vn{x), a; e M, n G N. Then iJ„ converges to H in norm resolvent sense 
as n — > 00, that is, 

Jiiin II (if„ - -2/L2(R.rf^))"' -{H- -2/L2(R.rf^))"'|| =0, zG C\M. (3.72) 
This follows from the resolvent identities, 

z G C\a{Hn), (3.73) 

- {H^°^ - ZlL2{R;dx))~^v[lL2{R;dx) " K {z)]'^ U {H^^^ - Z I {R;dx)) ' ^ , 

z G C\a{H), (3.74) 

where 

Kr,{z) = UniH^"^ - zlL2(^,dx))~^Vn, z G CV(ifW), n G N, (3.75) 
and K{z) is given by (3.26), and the fact that (3.71) implies 

Jim ||(«„ - ^;)(ifW - zh2^^.,dx))-"\^^^.^^.^,^^^ =0, ze CV(fW) (3.76) 

(cf. the detailed discussion in [13]). Thus, the spectrum of Hn converges to that of 
H as n ^ (X. In particular, for n G N sufficiently large, Hn has a simple eigenvalue 
Xn in a small neighborhood of z = satisfying 

Xn 0. (3.77) 



n— >oo 



We denote by 'i'n the corresponding eigenfunction of iJ„, associated with the eigen- 
value Xn of Hn, Hn'^n — A„*„, n G N. is then a constant multiple of the 
solutions satisfying (3.43) with z = A„ and V replaced by Vn- We may choose the 
constant multiple in such that 

lim ||*„ - *o|U2(M;da;) = 0. (3.78) 

In addition, we recall that (3.43) also implies that '^n and are exponentially 
bounded in x G M with bounds uniform with respect to n G N. 
In addition, we abbreviate 

= $„ = CnSn^n, ^n) L2(M.;dx) = 1' (^•''^) 

P„ =($„,•)$„, (3.80) 

Snf = sgn(W„)/, / G L2(R; dx), (3.81) 
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/ dxWn{x)\^n{x) 



(3.82) 

Kn{z) = Ko,n + Ki,niz-Xn) + O{{z-Xnf), (3.83) 

n G N, (3.84) 

and recall that 

- Un{H^°^ - \nlL'^{tL-dx)) ^«n*n = ^n, -Vn{H^°'' - XnlL^(R-dx)) ^U„$„ = 

(3.85) 

Hn^n = {H^°^ + Wn)<i^n = A„*„, -{H^°^ - A„7i2(M;d,)) " V„*„ = n G N. 

(3.86) 

We note that since m„, w„, Wn, n E N, arc all bounded operators on L^(R; dx), no 
operator closure symbols in Un[H^'^^ — XnlL^(]SL-dx)) ^^nj k = 1,2, are needed in the 
following computation leading up to (3.87). 
Next, one computes 

<iGip„L^{R;dx){PnKi^nPn) = -detp^i2 (k.^j.) (P„M„ (i7 - XnlL^(R;dx)) ^ VnPn) 

= - A„/l2(r.^^)) ^U„$„, - Xnl L-HVi;dx)) ^"«'^'n)L2(R;da;) 

= -Cn{{H^'°^ - XnlL^(R-dx)) Wn^n, {H'-°^ - XnlL2(R-dx)) ^n*n)i2(R.rf^) 



= -c„ 



lL2(M;da;)' 



n G N. 



Since 



Hm C„ = Co = 



/ dxW{x)\^o{x)(" 
Js. 



(3.87) 



(3.88) 



one finally obtains, 



detpoi,2(R.(j3.)(PoifiPo) = -detpgL2(s^.dx) 

\ J L^(R;dx) 

= - lim {^ri,Un{H'-°^ - XnlL'^{K-dx))~'^Vn^n] 



-C'o||*o||i2(R.d^) 

- / dxW{x)\^Q{x)\'^ 

-\\^0\\l2^R.dx)/i'"^0,U-^0)L^{R;dx) 



ll*0||i2(R.d^) 



|V1.M,2 



+ v;o||*ol 



ll*o|| 



(3.89) 

(3.90) 

(3.91) 

> 0. 
(3.92) 
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Here we applied the quadratic form equality 

< ll*olli2(R.rf,) + K,c||*o|li2(R.rf^) = -(^;4'0,U«'0)L^(R;dc.) 

= - [ dxW{x)\^o{x)\^ ^^'^^^ 
Js. 

to (3.90), to arrive at (3.92). □ 

3.2.2. The first factor: A posteriori computation. Before we proceed to a direct 
approach to compute the first factor on the right-hand side of (3.64), 

deti,2(R.dj;)(7£,2(R.dj;) - Kq- Po), (3.94) 

we will next determine det£,2(R.(j3.)(7j^2(R.(j3.) — Kq — Pq) by using the final answer 
(3.51) for ^■•(0). 

Theorem 3.7. Assume Hypothesis 3.2 and suppose € ad{H). Then, 

<ie^L2{M;dx){lL2(JR;dx) - Pq - Kq) 

= ^m [ dx[V{x)-Voo]^+{0,x)^-{0,x) (3.95) 
2Voo J«- 

= -[2yVr'll*±(0)|IZW)Ol*±(0)ll'=(K;ax) +K.||*±(0)||i2(R^,,)] 
x(*_(0),'J'+(0))i2(R.rf,). (3.96) 

{Here the equations for the + and — sign should be read separately.) 

Proof Combining (3.44), (3.51) (setting = 0), (3.64), (3.68), and (3.89), and 
taking into account that for some constants C± G M.\{0}, 

^o{x) = C±*±(0, x), xgR, (3.97) 
then yields for the first factor in (3.64), 

J^*(0) 

detL2(R-dx){lLHR;dx) - Pa- Ko) = — (p p\ 

= Co-^||*oir'^I7^ / dx'^.{0,x')^+{0,x') 

= [2^0^/'] II *0 II Z2(R;d.) (t'^'O, U*0)L2(K;d,) (*- (0), ^'+(0))i2(R.d,) 

= ^T72 I dx[V{z)-VJ^^+{Q,x)^_{Q,x) 

2V(X) 

= -K/^]"'ll*±(0)IIZW)[ll*±Wlli=(M;dx) + K.||*±(0)||i2(R,,,)] 

x(*_(0),'J'+(0))i2(R.rf,). (3.98) 

□ 

3.2.3. The first factor: Direct computation. Next, we proceed to a direct approach 

to compute the first factor on the right-hand side of (3.64), det2,2(ji.(;j.)(//^2(]5.^3.-) — 
Kq — Pq). This will now be an ab initio calculation entirely independent of the 
result (3.61). 
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Theorem 3.8. Assume Hypothesis 3.2 and suppose S ad{H). Then. 
deti2(K.rf^)(7i2(K.d^)-Po-i^o) = ^i2 / dx[V{x)-V^]e^^^'"''i>±{x), (3.99) 
where i/'± are defined by 
^±{x) = -*±(0,x) - ^ dx' sinh (yV2(2. _ x'))[V{x') - Vo.]i>±{x'), 

xeR. (3.100) 

{Here the equations for the + and — sign should be read separately.) 

Proof. Our strategy is to apply formulas (3.9) and (3.12) in [13] to the Frcdholm 
determinant detL2(jg^.^^^'^{Ij^2(^^.^^'^ — Pq — Kq) by appealing to the semi-separable 
nature of Pq + Kq upon an appropriate reformulation involving block operator 
matrices. To this end we introduce 



/l(^) = 


(-M(a;)e-^~'^ u{x)^o{xyj , 


(3.101) 


9i{x) = 


([2(Ko)i/2]-i«(x)e^i^'- Covix)^oix)y . 


(3.102) 


f2{x) = 


(-u(a;)e^~'^ u{x)'^o{x)) , 


(3.103) 


92{x) = 


([2(Ko)'/']-'^^(ar)e-^i''- Cov{x)'^o{x)y 


(3.104) 



and note that Pq + Kq is an integral operator with semi-separable integral kernel 

ru , ^/ /^ j h{x)gi{x'), x'<X, fo,nK\ 
{Po+Ko){x,x) = < (3.105) 
\f2{x)g2{x'), x' > x. 



In addition, we introduce the integral kernel 



sinh(li'-(.r-3;')) 



H{x,x') = fi{x)gi{x') - f2{x)g2{x') = u{x) ^ ' " v{x'), (3.106) 

and, temporarily assuming that 

supp iy — Vao) is compact, (3.107) 
the pair of Volterra integral equations 

/•oo 

/i(x) = /i(x) - / dx' H{x,x')h{x'), (3.108) 



f2{x) = f2{x) + f dx' H{x, x')f2{x') (3.109) 

J — OO 

for a.e. a; e M. Applying Theorem 3.2 of [13] (especially, (3.9) and (3.12) in [13]) 
one then infers 

deti,2(R.rf3.)(7i2(R.rf3.) - Po - Kq) = detc2 - j dx gi{x)fi{x)j (3.110) 

= detc2 (^l2 - J^dxg2{x)f2{x)j , (3.111) 



22 F. GESZTESY, YU. LATUSHKIN, AND K. ZUMBRUN 

with I2 the identity matrix in C^. Introducing \l/± as the solutions of the pair of 
Volterra integral equations 

- ^ / sinh {V^\x - x')) [V{x') - V^\^±{x'), a; e M, 

(3.112) 

a comparison with (3.108) and (3.109) yields 

fi{x) = -u{x)^+{x), f2{x) = -u{x)^-{x) for a.e. x eR. (3.113) 

Writing 

^±{x) = {^±{0,x) ip±{x)), xgR, (3.114) 

(3.112) yields {x G M) 

i>±{x) = -^oix)-^ / dx' smh{V^/\x-x'))[V{x')-V^]i>±{x'), (3.115) 
and 

-1 /-iboo 

*±(0,ar) = e^^-'- - ^ / '^a;' sinh (yV2(^ _ x'))[V{x') - V^]^±{0,x% 

(3.116) 

in accordance with (3.43) for 2 = 0. Because of (3.97), ^q{x) = C±^'±(0,a;), a; € R, 
one infers that 

i)±{x)=C^i)±{x), a;eM, (3.117) 

with satisfying (3.100). 
Next, one computes 



, oxp(vi/^x) cxp(vi/^.) 



d2{x)h{x) 



-Co[V(x)-Voo]*o(x)*+(0,x) -Co[V(x)-Voo]*o(a:)j/;+(a:) 

^xp(-V-4''^.), 



-Co[V(^)-Vool*o(2:)*-(0,a;) -Co[y(a;)-ycx,l*o(a;)'0- (x) 



(3.119) 



Using the fact that G ad{H), and hence deti2(R.<j3.)(7j,2(R.rfj.) — ii'o) = 0, one 
obtains from taking 2: = in (3.47) and (3.50), 



detL2(M.,dx)(/L=(M;dx) -Ko) = l + ( e±^~ [F(a;) - Foo]*±(0,a;) = 0. 

(3.120) 



21/00' -'R 

Thus, (3.110) and (3.111) together with (3.118)-(3.120) yield 



deti2(K.d^)(7i2(K.d^) - Po - i^o) = -^oj^dx [V{x) - yoo]*o(a;)*±(0,a;) (3.121) 

1^ / dxe±^i''-[l/(x)-14,]V'±(x) 



X 

2Vi 
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(where the equations for the + and — sign should be read separately). Applying 
(3.39) for Xj = and (3.97) one obtains 

detL2(M;dx){lL^(R;dx) - Po - Kq) = — — / rfx e=^^~ [T/(a;) - Voo]^P±ix) 

2Voo Jr 

(3.122) 

using (3.117) in the last line. 

To remove the temporary compact support assumption (3.107) we first note that 
by a standard iteration argument, the Volterra equations 

e±^i''-V±W = -e±^~'-vl/±(0,x) T / dx' [e^^y'J'i--') - l] 

IVoo Jx 

X [Vix') - Kc.]e±^i'""'V^±(-^'), xGR, (3.123) 
have unique and bounded solutions on M, which together with their first derivatives 
are locally absolutely continuous on M, as long as the condition (3.16) is satisfied. 
This follows since there exists a constant C > such that [e^^-^ ^^^'±(0, a;)| < C, 
X G M. Thus, the right-hand side of (3.99) remains well-defined under condition 
(3.16) on V. 

Next, similarly to the proof of Theorem 3.6, we choose compactly supported 

sequences Un,Vn G L^(]R: dx), n £ N, such that Wn = Vn — Voo = UnVn converges 
to W = V — Voo = uv in L^{R;dx) as n ^ oo and introduce the maximally 
defined operator i?„ in L'^{R;dx) associated with the differential expression i„ = 
— {d'^/dx'^) + Vn(x), X G M. Since _ff„ converges to H in norm resolvent sense 
(this follows in exactly the same manner as discussed in the proof of Theorem 3.6), 
the spectrum of iJ„ converges to that of as n ^ oo. In particular, for n G N 
sufficiently large, _ff„ has a simple eigenvalue A„ in a small neighborhood of z = 
such that A„ — > . Multiplying Vn by a suitable coupling constant g„ G M, 

n — >oo 

where gf„ — > 1, then guarantees that the maximally defined operator iJ„(g„) in 

n— ^oo 

L^{R;dx) associated with the differential expression ^^(^n) = —{d!^/dx^)+gnVn{x), 
a: G M, has a simple eigenvalue at z 0, in particular, G (yd{Hn(gn))- (Multiplying 
V by Qn changes the essential spectrum of Hn{gn) into [gnVoo, oo), but since — > 
and 5f„ — > 1 as n — > 00, this shift in the essential spectrum is irrelevant in this proof 
as long as n G N is sufficiently large.) 

Finally, the approximation arguments described in the proof of Theorem 4.3 of 
[13] permit one to pass to the limit n — > oo establishing (3.99) without the extra 
hypothesis (3.107). □ 

It remains to show that the two results (3.95) and (3.99) for the Fredholm de- 
terminant dct^2(-jj.^2;)(-^L2(R;d2:) ~ -Po ^ Kq) coincidc. This will be imdertakcn next. 

Theorem 3.9. Assume Hypothesis 3.2 and suppose G ad{H). Then the expres- 
sions (3.95) and (3.99) for deti2(R.rf^)(7j^2(R.rf^) - Po - ^o) coincide. 

Proof. To keep the arguments as short as possible, we first prove that (3.95) and 
(3.99) coincide under the simplifying compact support assumption (3.107) onV — 
Foo- Again, the general case where V satisfies Hypothesis 3.2 then follows from an 
approximation argument. 
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More precisely, we suppose that 

supp {V -Voo) CI [-R, R] for some i? > 0. (3.124) 
Proving that (3.95) and (3.99) coincide is then equivalent to showing that 

/R pR 
dx [V{x) - Ko]e=^^~ '^V±(a;) = dx [V{x) - Voo]*+(0,a;)*_(0,a;). 
-R J-R 

(3.125) 

We start with the right-hand side of (3.125): First we note that 
^Wr(V'±,*q.(0))(x) = [V{x) - V^]^+{0,x)^_{0,x) for a.e. x eR, (3.126) 

where we used that 

- V±(a;) + V{x)tp±{x) = [V{x) - V^]'^±{0, x) for a.e. x€R, (3.127) 
which in turn follows by twice differentiating (3.100), and 

- *±(0,a;) + V{x)^±{0,x) = for a.e. xgR. (3.128) 
Thus, one concludes that 

dx[V{x)-Voo]-^+{0,x)-^-{0,x) 

R 

d 

= / dx—WT:(ip+,^(0))(x) 
J-R dx 

= Wr(V;+, * -(0))(ii) - Wr(V+, *-(0))(-i?) 
= -Wr(f+(0),*_(0))(i?) - Wr(^+,*_(0))(-i?) 
= -Wr(V'+,*-(0))(-i?). (3.129) 
Here we employed that 

tp+{x) = -*+(0,a;) for a; > ii (3.130) 

(cf. (3.100)) and 

2V^^J='{0) = Wr(*+(G),*_(G)) = (3.131) 
since by hypothesis, £ ad{H). Similarly, one computes 

i-R 

dx[V{x)-Voc]^+{0,x)^-{0,x) 

R 

d 

= / da;— Wi-(?,^'_.*+(0))(a;) 
J-R dx 

= Wr(V'-, *+(0))(i?) - Wr(7/'_, *+(0))(-i?) 

= Wr(V'- , *+ (0) ) (i?) + Wr (*_ (0) , *+ (0)) (-i?) 

= Wr(V;_,vl/+(0))(i?), (3.132) 

where we used 

ijj-ix) = -^--(0,2;) for x<-R (3.133) 
(cf. (3.100)) and again (3.131). In particular, one concludes that 

Wr(V'-,*+(0))(i?) = -Wr{ip+,^-{0)){-R). (3.134) 
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To compute the left-hand side of (3.125) we first note tliat 

= -[V{x) - Voo]e=^^~ '"^Vila;) for a.e. a; G M, 
where we employed again (3.127) and (3.128). Thus, one infers that 
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da; [T/(x) - Ko]e^- "V+(^) 



R 



dx 



R 



d_ 

dx 



Wr(v+(a;),e^~'") + ^Wr(«'+(0, x), 



: Wr 
+ 



(V'+(a;),e^~ 



Wr 



x=R 



(^+(x),e^i^^- 



x=-R 



■ Wr(*+(0, a;), | _^ - Wr(*+(0, x), 

= -Wrf^'+(0,x),e^~'^) + Wr(l'+(0,x),e^'i 

V / x—R V 

^+{0,x),e^^''^) 



x=-R 



x=R 



Wr 



x=-R 



- Wr(v+(a;),e^~'""'^ 

= -Wr(?A+.*_(0))(-i?) - Wr(4'+(0),*_(0))(-i?) 
= -Wr(VJ+,*-(0))(-i?). 

Here we used again (3.130) and (3.131) as well as 

^v^''x ^ ^_^o,a:) for x < -R 

(cf. (3.116)). Similarly one computes 



x=-R 



R 



dx [V{x) - Foole-^- '"^-(x) = Wr(V'-, *+(0))(i?). 



Taking into account (3.134), this completes the proof of (3.125). 



(3.135) 



(3.136) 
(3.137) 

(3.138) 
□ 



3.3. A formula of Simon. Finally, we turn to an interesting formula for the Jost 
solutions ^'±(z, • ) in terms of Fredholm determinants derived by Simon [37]. 

To set the stage, we abbreviate M± = (0, ±oo) and introduce the one-dimensional 
Dirichlet and Neumann Laplacians perturbed by the constant potential V^o, H^^'^ 
±;dx) by 



and in 



Ff'^ = -A_ + v^^, dom(fff'^)=Fo2(R±), 



(3.139) 



K 



(0),N _ 



± + dom{H^^>''')={geH\R±)\g'{0±)=0}. (3.140) 

Next, we recall that 

^±{z,0) = detL2(M^:dx){l + uiH^±^''' - z)-'vy lm{iz-Voo)'^^)>0, (3.141) 

a celebrated formula by Jost and Pais [20] (in the case Voc = 0). For more de- 
tails and background on (3.141) we refer to [13] and the references cited therein. 
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Moreover, it is known (cf. [12], [13]) that 

0) = ±i{z - V^f'^deiL^K±,dx) (l + '"^ - z)- .J , ^^^^ 

lm{{z-V^f'^) >0. 

We conclude this section by presenting a quick proof of the representation of the 
Jost solutions 'i!±{z,x) and their a;-derivatives, ^'j^{z,x), in terms of symmetrized 
perturbation determinants, starting from the Jost and Pais formula (3.141) and its 

analog (3.142) for *±(z,0): 

Lemma 3.10 ([37]). Suppose V satisfies (3.16) and let Im((z - K>o)^^^) > 0, 
a; e R. Then, 

*±(z,:r) = e±'(^-^~)'''- 

X detL2(R^.rf3,)^/i2(R^.rf3,)) +'u(- +x)(ij£°^'"° - z)"^w(- +x)^, 

(3.143) 

^'±{z,x) = ±i(0- yoo)^/^e±'(^-^~)'''^ 

xdetL2(^^^.dx)(^lL^{K±;dx))+u{ - + x)(ij£°^'^ - 2;)"^?;( • +x)y 

(3.144) 

Proof. Denoting Vy{x) = V{x + y), x,y G M, and by \l/y_±(z, •) the Jost solu- 
tions associated with Vy, an elementary change of variables in the Volterra integral 
equation (3.43) for 5'j,,± yields 

^y,±iz,x) = eT^(^-^-)^''^*±(z,x + y), 

1/2 (3.145) 
*;,+ (z.,7:) =eT^(^-^~) y^'+{z,x + y). 

Taking x = in (3.145) implies 

*^,±(^,0) = eT*(^-^-)'''^*±(^,t/), (3.146) 

*;,±(^,0) = e^*(^-^~)'''^*'±(^,y). (3.147) 

Using the Jost-Pais-type formulas 

^y,±{z, 0) = deti2(R^.rf,) (/ + «(•+ 2/)(iff '"^ - ^)-^?;( ■+?/)), (3.148) 

%^±{z, 0) = ±i(^ - Ko)'/'deti2(M^.rf,) (/ + u( • + t/)(iff '"^ - ^)-it;( • + y)) , 

(3.149) 

an insertion of (3.148) into the left-hand side of (3.146) proves (3.143). Similarly, 
an insertion of (3.149) into the left-hand side of (3.147) yields (3.144). □ 

4. The multi-dimensional case 

In the previous section, we have illustrated within the simple setting of one- 
dimensional scalar reaction-diffusion equations how the stability index may be 
equally well calculated from an Jost/Evans function point of view, or else, us- 
ing semi-separability of the integral kernels of Birman Schwinger-type operators, 
directly from first principles using Fredholm determinants. We conclude by describ- 
ing, again within the reaction-diffusion setting, an algorithm for multi-dimensional 
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computations via Fredholm determinants, based on semi-separability of the integral 
kernels combined with Galerkin approximations. 

4.1. Flow in an infinite cylinder. Consider a scalar reaction-diffusion equation 

wt = Aw + f{w), (4.1) 

on an infinite cylinder x = (xi,X2, ■ ■ ■ ,Xd) £ M x O, where A = is the Lapla- 
cian in the x- variables, f2 C is a bounded domain, w and / are real- valued 

functions, 

/ e C^{R). (4.2) 
In what follows we will assume that = [0, 27r]'^~^ and consider only the physical 
cases d = 2,3. Unless explicitly stated otherwise, we will always assume that 
periodic boundary conditions are used on the boundary dCl of fl (viewing O as a 
(d — l)-dimensional torus in the following) if d = 2, 3. For a; € M x J7 we will always 
write X = {xi,y), where xi G R and y = {x2, ■ ■ ■ ,Xii) € fl, and similarly, x' = 
{x'l,. . . , x'j) = {x[,y'), y' = (x2, . . . , x'j). We will abbreviate dx = dxidx^ ■ ■ ■ dxd 
and dy = dx2 ■ ■ ■ dxd, and frequently use the fact that the space L^(IR x 17; dx) = 
L'^(R;dxi;L^{Q;dy)) = L'^(Q;dy;L^{M.;dxi)) is isometrically isomorphic to the 
space L'^{R; dxi)) via the discrete Fourier transform in the y-variables: 

w{x)= Wj{xi)e'^-y, x={xi,y)GRxQ, (4.3) 

where 

Wj{xi) = {2ny-'^ I dyw{xi,y)e-'^ y, Xi^R. (4.4) 
Jn 

4.1.1. Galerkin-based Evans function. We first review the; Galerkin approach de- 
scribed in [23], in which a standard Evans function is defined for a one-dimensional 
truncation of the linearized operator about a standing-wave solution in a series of 
remarks. 

Remark 4.1. Under our standing assumption of periodic boundary conditions on 
dQ, there exist planar steady-state solutions U = U{xi), where U is the solution 
of the corresponding one-dimensional problem (3.4) described in Section 3. Lin- 
earizing about U, cf. (3.1)-(3.9), denoting V{xi) = —f'{U{xi)), and taking the 
discrete Fourier transform in directions X2, ■ ■ ■ ,Xd, one obtains a decoupled family 
of one-dimensional eigenvalue problems 

= {Lj - A)^ = - bf - A - V{x^)y, (4.5) 

indexed by Fourier frequencies j = (72, . . . ,jd) G Z''^^, each of which possess a 
well-defined Evans function and stability index. At j = and A = 0, there is an 
eigenfunction U'{xi) associated with translation invariance in the xi-direction of 
the underlying equations; for other j, there is typically no eigenfunction at A = 0. 
Asymptotic analysis as in [1], [30] yields a trivial, positive stability index for jjj 
sufficiently large, so that computations may be truncated at a finite value of 

Remark 4.2. In the above example, the operators Lj are real- valued (i.e., map 

real- valued functions into real- valued ones), hence a stability index makes sense. 
For more general, non-self adjoint operators, one may expand in sines and cosines 
to obtain a family of real- valued eigenvalue equations for which a stability index may 
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again be defined. This principle extends further to generahzed Fourier expansions 
in the case of general fi, requiring only real-valuedness (in the above sense) of the 
original (multi-dimensional) operator L. 

Remark 4.3. More generally, consider a standing-wave solution U = U{x) that is 
not planar, but only converges as xi ±00 to a constant state Uoo- We assume 

(xi ^ ||C/(xi, •)llff3/2(n)) G (iinL°°)(R;da;i). (4.6) 

This writing means that the map Xi ^ U{xi, ■ ) from R into the fractional Sobolev 
space H^^/'^{Q) is both an L^- and L°°-function with respect to the variable xi. 
Linearizing about U, denoting by = ■^(a;) the corresponding eigenfunction, and 
taking the Fourier transform in directions a;2, . . . , Xd, one obtains a coupled family 
of one- dimensional eigenvalue problems 

o=(^-bf -A^v^i-Cn^i, •)*^(.))(i). iez'^-i. (47) 

Here * denotes convolution in j , 

iIj{x)= iPj{xi)e'^-y, x = {xi,y) €Rxn, (4.8) 

and V{xi,j) denotes the value of the Fourier transform of V{xi, ■ ) in the variable 
y = {x2, ■ ■ ■ ,Xd). Following the approach of [23], one may proceed by Galerkin 
approximation, truncating the system at some sufficiently high-order mode \j \ ^ J, 
to obtain again a very large, but finite, real- valued eigenvalue ODE in xi , for which 
one may define in the usual way an Evans function and a stability index. 

Remark 4.4. While we focus primarily on periodic boundary conditions on Oft 
throughout this section, one can treat other boundary conditions such as Dirich- 
let, Neumann, or more generally, Robin-type boundary conditions in an analogous 
fashion. The key fact used in (4.3) and (4.4) is the eigenfunction expansion associ- 
ated with the discrete eigenvalue problem of the self-adjoint Laplacian in L^{fl; dy)) 
with periodic boundary conditions on d^. The latter can be replaced by analogous 
discrete eigenvalue problems of the Laplacian with other self-adjoint boundary con- 
ditions on dVl. 

4.1.2. Fredholm determinant version. We now describe an alternative method based 
on the Fredholm determinant, in which the Jost and Evans functions are prescribed 
canonically as characteristic determinants, but computed by Galerkin approxima- 
tion: that is, we approximate the determinant rather than the system of equations. 

Specifically, consider again the general situation of Remark 4.3 of a solution U 
of (4.1) decaying as xi ±00 to some constant state ?7oo- As in (3.9), we denote 
V{x) = —f'{U{x)), a; G ri X R, and assume that 

Voo = -f'iUoo) > 0. (4.9) 
We define the Birman-Schwinger operator K{z) similarly to (3.26), (3.22) as 

K{z) = -u{hIII - zlL2(K^n;d.))''v, z € C\[T4o, 00), (4.10) 

with iJo'p the self-adjoint realization of the differential expression — A^, -|- V^o in 
L'^{R X Q;dx) with periodic boundary conditions on dfl, and 

u{x)=sgn{V{x)-V^)v{x), v{x) = \V{x)-Voo\^/^ (4.11) 
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for a.e. x £ M x fi. 

For later reference, we define also the asymmetric rearrangement of K{z) by the 
formula 

IC{z) = -{H^I^- zlL2(M.xn;dx)) ^uv 

= -{HSi-zlL^{Rxa;d.)y\v-V^), zeC\[V^,oo). (4.12) 

Then, under the assumption 

{xi ^ \\V{xu ■ ) - Voo \\Looin;dv)) e {L^ n (M; dx^), (4.13) 

we have the following result generalizing the one-dimensional case [11, Lemma 
2.9]. Fix z e C\[Foo,oo). Passing to adjoint operators, if needed, with no loss of 
generality we will assume below that Im(2;) ^ and fix the branch of the square 
root such that Rc((t4o + |jP ~ -z)^^^) > for each j e a choice consistent 

with the choice of Q in (4.95). 

Lemma 4.5. Assume (4.13). // dimensions d ^ 3 then K{z), IC{z) G S2(i^(R x 
il;dx)) for each z e C\[14o,oo). Moreover, the condition on the dimensions is 
sharp. 

Proof. The operator Hq^^, since constant-coefficient, decouples under the Fourier 
transform in the variables X2, ■ . ■ ,Xd. Consequently, the integral kernel of the resol- 
vent of -ffo°p' denoted by {H^^^ - zlL2(^Rxa;dx)) ^{x,x'), may be found exphcitly 
by a Fourier expansion and, using (3.27), can be expressed as a countable sum of 
scalar integral kernels: 



= J2 i^-Voo- |j|2)-i/2ei(-^oo-bf )^'^|x.-x',|gii.(^-^'), (4.14) 

X = {xi,y),x' = {x[,y') eRx fl, 

where y = {x2,. . ■ , Xd) y' = {x'2, . , . , x'^) £ fl, and j € Z''"^ denote the Fourier 
wave numbers in these directions. Using Parseval's identity, we obtain for any fixed 
a;' e M X O that 



t2(RxQ;dx) 

2 



= i E \z-V^-\j\X"\s\n\{^rg{z-V^-\jf))y\ (4.15) 



where Im((z - - |jp)^/^) > for 2 e C\[Ko,oo) due to Im(z) ^ 0. Since 
arg [z — Voo — IjP) ^ tt as |j| ^ 00, there is a constant c = 0(2;) such that 

||«) -./^.(Mxn;.x))"\-,a:')||l.(«,n;.x)^c E (^oo + b•^-^/^ (4.16) 
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hence is finite if and only if d ^ 3. We recall the formula for the Hilbert-Schmidt 
norm of the Hilbert Schmidt operator K with integral kernel K(x,x') (see, e.g., 
[38, Thm. 2.11], [46, Sect. 1.6.5]): 

\\K\\B2{L^(V.xn;dx)) = , • )IU2((Kxa)x(Rxn);da;dx')- (4-17) 

Using (4.17) and (4.16) to estimate the integral kernels of (4.10) and (4.12), one 

infers 



\K{z,-. 



L2((Kxn)x(Mxa);dxdx') 



< ||«(-)||«!.-^^i^(Mxa;dx))"'(-, •)||^.(Mxa;dx') ^.(KxQ;d.) 



lL'»(Kxa;dx') 



<c (y^ + |j|2)-3/2 . (4.18) 



((Rxn)x(Rxn);dxdx') 
1 1 2 



L2(MxO;da;') 



n,p ^-^L^(M.xn;dx)) \ , f\\L^(M.xn;dx)' 

<c Y (Kx> + lil')-'/' • ||wt^||i.(Mxa;dx')' (4-19) 



and finds that K{z) and /C(z) are Hilbert-Schmidt operators for d ^ 3, as claimed. 
In the decoupled case, where u = u{xi), v = v{xi), these estimates are sharp, 
showing that in general K{z), IC{z) are Hilbert-Schmidt only for d ^ 3. □ 

Definition 4.6. Assume (4.13). Generalizing the one- dimensional case (3.50), we 
introduce in dimensions d= 2,3, a 2-modified Jost function defined by 

J^2{z) = det2,mRxU;dx){lL^Rxn;dx) " z G C\[Voo, oo). (4.20) 

By the determinant property (2.42) one then obtains 

J^2{z) = det2,L'^(Rxn;dx){lL^(Rxn;dx) " K{z)), z e C\[Vao, oo), (4.21) 
which could equivalently have been used to define ^2{z)- 

Remark 4.7. Similar to Lemma 4.5 calculations show that K{z) belongs to suc- 
cessively weaker trace ideal classes as d increases, for which a higher-modified Jost 

function may be defined as a higher-modified Fredholm determinant. We restrict 

our attention here to the main physical cases d = 2, 3. 

Remark 4.8. Comparing with formula (3.50) for the non-modified Jost function 
given in the one- dimensional d = 1 case (when K{z) and )C{z) are trace-class 
operators), (2.41) implies the relation 

det2 i,2(R^^ )(/j^2(R ^;^ ) - K{z)) 

(4 22) 

= deti2(K.rf,,)(7^.(K.rf,,) -i^(^))e*'--^*<^»i)(^(^». 

Thus, (4.20) differs from (3.50) by a nonvanishing analytic factor e*''^'(«i<i-i)^^^'^'\ 
and hence for practical purposes the use of dct2,L^(^-dxi){ ' ) dcti2(j^.d,j.^^( ■ ) in 
(4.20) for d = 1 are equivalent. For d > 1 we make a related comment in Remark 
4.14. 



DERIVATIVES OF (MODIFIED) FREDHOLM DETERMINANTS 



31 



4.1.3. Galerkin approximations. Next, we approximate J-2 by a Galerkin approx- 
imation, working for convenience with the asymmetric version (4.20). We will 
augment (4.13) with the more restrictive, but still typically satisfied, condition 

(xi ^ \\V{xu •) -^oo||H3/2(n)) e L\R;dXi). (4.23) 

As in (3.17), we introduce 

W{x) = V{x) - Kc- xeRxn, (4.24) 
and expand W into a Fourier scries in variables y = (x2, . . . ,Xd) so that 

W{x)= Wm{xi)e"" y, x = {xi,y) eRxn. (4.25) 

Substituting (4.14) into (4.12), we obtain an expansion of the integral kernel of the 
operator )C{z): 

IC{z,x,x') = - ^^^e^^-y {V{x') - V^)e-^= ' y (4.26) 

e-<~y^+\3\^-zY'^\x^-x'^\ _ 

Introducing 

fi{x) = -2-^e'i-y{Voo + bf - ^)-V2e(V'oo+bf -.)^/=xi^ 
fiix) = -2-^e'i-y{Voo + bf - ^)-i/2e-(Vo„+bf -.)V=xi^ 
g{{x') = e-(^~+ljl'-^)'''^'i(F(a;i) - V^)e-'^-y', 
gi{x') = e(^~+l^l'-^)'''^'i(F(a;'i) - Foo)e-'^' 
we obtain from (4.26) an expansion of lC{z, a countable sum 



(4.28) 



IC{z,x,x') = /^J6Z<i-i f({x)gi{x'), x[ > xi, 
[Jljez"-^ f2ix)9i{x'), xi>x[, 



(4.29) 



of scalar integral kernels that are semi-separable in xi. 

Truncating (4.27) at some finite wave number J or, equivalently, Fourier expand- 
ing ff., g^. in (4.29) in variables y and y' and truncating the resulting series at some 
finite wave number J, we obtain a sequence of Galerkin approximations 

)Cj{z,x,x') = - ^|,|2 ,u/2 W-,-m(xDe— ^ (4.30) 

E e^J -^(/?)^.(xi)(,gl)_^(xl)e-'™-''', x\> x^ 

_ j ^ _ (4 31) 

E e^^-^(/ll(^i)(5^)_„W)e-™-^', a;i>a;i, 

where (fl)^ denote the Fourier coefficients of 

m^uy)= E ifi)^e^"'-' = {7l)/'-', (4.32) 
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and 

A)J^i) - %™K)e(-^)'="^(^~+l^1^-)^^^-S fc = 1,2, (4.33) 
denote the Fourier coefficients of the function 

gi{x') = e^-i^'^'^^^+l^'I'-^^'^'^'^WCajOe-'^-"', k = l,2. (4.34) 

We denote by lCj{z) the integral operator on L^(R x SI; dx) with the integral kernel 
(4.30), (4.31). 

Theorem 4.9. Let z G C\[Foo,oo). Then under assumptions (4.13) and (4.23), 
JCj{z) e S2(L^(R X fi; dx)) for dimensions d = 2, 3. Moreover, for d = 2,3, ICj{z) 
converges in the Hilbert-Schmidt norm to JC{z) at rate j('^~^y^ as J oo and 
hence the sequence J^2,j{z) defined by 

^2,j{z) = det2,I,2(Rxa;dx)(-fL2(RxQ;da;) " ^j{z)), (4.35) 

converges to J^2{z) = det2,i,2(Rxn;dx)(-'"i,2(Mxa;dx)-^(-2)) as J ^ oo at rate J^'''-^^'^ . 
Proof. Using (4.27) and (4.30), one obtains 

J oc 



lC{z,x,x')-ICj{z,x,x') = -(j2 E + E E) 



(4.36) 



j|=0|m| = J+l lj| = J+l |m|=0 ' 



By the triangle inequality and Parseval's identity one therefore infers 

-\\)C{z) - ^j(^)||i2((Kxn)x(Mxa);dxxda:') 



J OO 



E E (•) 

|j|=0|m|=J+l 



OO OO 



L2((Rxn)x(Mxa);dxda;') 
2 



E E(-) 

|j|=J+l |m|=0 

J CxD 

E E 

|j|=0|m|=J+l 



i2((Mxa)x(MxQ);dxdx') 

-l./l"-')' '-:-'-i-2:;L 



OO OO 

E E 

|j|=J+l |m|=0 



2(K,o + liP - ^)V2 



W^.-m(a^'i) 



I,2(RxM;dxidx'i) 
2 



2(Ko + |iP-^)V2 



(4.37) 



(4.38) 



I,2(RxR;dxidx'i) 



We will now estimate (4.37) and (4.38) separately. Using arguments similar to 
(4.15)-(4.16), one observes that the sum in (4.37) can be estimated as follows: 



(4.37) = J2 



b1=o 



2(^00 + |jf - ^)^/2 



E 11^^— (•)! 



<C Y: (Ko + Un-^/^ E ll^'»(-)|ll^(«;.xi) 
jeZ<*-i |m| = J+l 

|2 



(4.39) 
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In the last equality we used that the series X^^g^d-i converges due to 1 < c? < 3, 
and Parseval's identity for W — Wj, where 

Wjix)= Wmixi)e^"'-y, x={xi,y)€Mxn, (4.40) 

is the truncation of W. By the Sobolcv embedding W^^^'^iQ) ^ W^/'^-'^{n) (cf., 
e.g., [19, Theorem 1.6.1], [42, p. 328, Eq. (8)], [44, Sects. I.4-I.6]) and a standard 
Cauchy-Schwartz argument, one infers for each xi G R, 

\\w{x„.)-wj{xu-)\\i^n;dy)= E \^jix^)\'= E ui'/'l^.-(^i)l'br'^' 

\J\>J \J\>J 



< ( E mwM)\") ( E 



1/2 

3 



<c||W^(a;i,.)ll|3/.(o)-/^'-'^/', (4.41) 

Here we used standard notation M^*'^(-) = H.^{-) for Sobolev spaces. Thus, by 
(4.39), 

(4.37) <t^il|W^lll.(,^,,^^,3.(,/^''-^^^^- (4-42) 

Likewise, similarly to (4.15)-(4.16), one observes that the sum in (4.38) can be 
estimated as follows: 



<c £ (v;, + |,n-3/2 g 

b1=J+l |m|=0 



2 

L2(R;da:i) 



|2 



<c' ^ Ur'l|W^lli^(Exn;d.') ^^2J''-^llW/|li.(K,,,^,,,). (4.43) 
b' 1=^+1 

Combining (4.37), (4.38), (4.42), (4.43), and using [46, Section 1.6.5] as in (4.17), 

one arrives at the estimate 

ll^-^j|li2((Mxa)x(Rxn);dxdx') ^ 2||W^|||2(jixQ;da;')(*^l-^^''~^^^^ + '^2^''"^), (4.44) 

yielding the claimed result. □ 

Next, we take a closer look at properties of the integral operator K,j{z) in i^(]R x 
il; da;) with integral kernel given by (4.30), assuming at first that 

T?,_™ e L''{^;dx{), m,j G Z, |m|, |i| ^ J. (4.45) 
Using the fact that i^(R x fl;dx) decomposes into 

L'^{R X Q; dx) = L^(M; dxi) O L'^{Q; dy), (4.46) 
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we will exploit the natural tensor product structure of the individual terms ICm,j{z) 
in 

JCj{z)= J2 '^rnjiz), zeC\[V^,oo), (4.47) 

\m\,\jKJ 

where lCm,j{z), m,j G Z'^~^, |m|, |j| < J, are integral operators in i^(R x Q,;dx) 
with integral kernels given by 

e-{Vo.+\j\^-zy'^\x,-x[\ 

z e C\[yoo,c^), m,j G Z^-\ |m|, |i| s$ J. 

With respect to the tensor product structure (4.46), the operators ICm,j{z) decom- 
pose as 

}Cm,jiz) = Amjiz) (g) Bm,j, Z £ C\[Foo, Oo), (4.49) 

where the operator 



^m,j{z^X,x) — e " of'T/ I '12 M/2 ^^j—mi^l)^ ^) 



•^m,j{z) = ( - (d'^/dxl) + {V^ + - ^;)/z,2(R.d^^)) Vj_„, z e C\[Foo, oo), 

(4.50) 

in L^(M; dxi) has the integral kernel 

An,j(^,a:i,a;i) = 2(Kx> + |j? - z)y^ Wj-m{x[), (4.51) 

and Bm,j in i^(0; dy) has the integral kernel 

Br„j{y,y') = e'^-ye-'^-y'. (4.52) 
In particular, each Bm,j is a rank-one and hence trace class operator on Z/^(0; dy), 

Bm,jGBi{L\n;dy)). (4.53) 
Next, Am,j{z) is a Hilbert-Schmidt operator on L^(M; dxi), 

Amjiz) eB2{L'^{R;dxi)), ^ G C\[Foo, oo), (4.54) 
if and only if (cf., [38, Thm. 2.11], [46, Sect. 1.6.5] and (4.17)) 

Wj-m€ L'^{^;dxi). (4.55) 

In addition, applying [38, Theorem 4.5, Lemma 4.7], Am,j{z) is a trace class oper- 
ator on (R; dxi ) , 

Am,j{z)eBi{L\R;dxi)), zeC\{V^,oo), (4.56) 

if and only if 

Wj-m e (i^(K; dxi)) . (4.57) 
Here the Birman-Solomyak space £^{L'^{m.;dxi)) is defined by 

1 /2 

e'{L\R;dx,)) = i^feLl,{R) E(_^ dx,\f{x,)f^ <(X)|, (4.58) 

with Qn the unit cube in M centered at n e Z. We recall that (cf. [38, Ch. 4]) 
L^{m; (1 + |a;i|)'^da;i) C f(L'^{R; dxi)) g L^{M- dxx) D L2(R; rfa;^) 
for all 5 > 1/2. ^"^'^^^ 
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We note in passing, that the symmetrized version Am,j{z) of Am,j{z), given by 

^m,j{z) = Uj-m{ - {<Pldx\) + (l/oo + \j? - z)lL^R;dxi))~^i'j-m, qq^ 

Z G C\[yoo,Oo), 

where 

Uj-mixi) = agn {W.i-,nixi))vj-rnixi), ^^-^(a;!) = | (xi ) | (4.61) 

for a.e. Xi G M, is a trace class operator under the weaker assumption 

Wj-m&L\R;dxi). (4.62) 

Given these preparations, we can now summarize Hilbert-Schmidt and trace 
class properties of )Cj{z) as follows: 

Lemma 4.10. Assume z ^ C\[Voo,oo). Then, 

{i) K,j{z) is a Hilbert-Schmidt operator on L^(M x f2; dx) ifWj-m S i^(M; dxi) for 

all m,j G Z, \m\, \ j\ ^ J. 

{ii) K,j{z) is a trace class operator on L^(R x Q; da;) ifWj-m G £^ [L"^ {R; dxi)) for 
all m,j G Z, |m|, \j\ < J. 

Proof. Since the sum in (4.47) is finite, it suffices to prove the Hilbert-Schmidt and 
trace class properties of ICmj{z) for fixed m,j. Since by (4.49), 

l^m J (Z) I = \^m,j (Z) I ® \Bm,j \ , (4.63) 

where as usual, \T\ = {T*Ty^'^, the singular values of K-mjiz) (i.e., the eigenvalues 
of \K.m,j{z)\) are square summable, respectively, summable, if and only if the sin- 
gular values of Amjiz) arc square summable, respectively, summable, since Bmj 
is a rank-one operator and hence has precisely one nonzero singular value. This 
follows from the well-known fact that the spectrum of a tensor product Ai ig) A2 
in Hi (8) H2 of bounded operators Aj in the complex separable Hilbert spaces Hj , 
j = 1,2, is given by the product of the individual spectra (cf., e.g., [32, Theorem 
XIII.34]), that is. 



(4.64) 



a{Ai 8)^2) =cr(Ai) • cr(^2) 

= {A G C I A = A1A2, X, G aiAj), J = 1, 2}. 

Thus, one can apply (4.54), (4.55), respectively, (4.56), (4.57). □ 



Of course, the condition on W in the Hilbert Schmidt context in Lemma 4.10 is 
much weaker than condition (4.23) since only finitely many Fourier coefficients Wk 
of W are in involved in the former, while the stronger condition (4.23) is used to 
prove the convergence of ICj in Theorem 4.9. 

Remark 4.11. More generally, any useful approximation of -ffj'j'p may be cm- 
ployed, not necessarily an eigenfunction expansion or one attached to a Fourier 
basis. In particular, in the case that Voo is not constant in X2, . . . ,Xd, one may 
proceed alternatively by Galerkin approximation to approximate Hq^^ as the limit 
of operators with semi-separable integral kernels corresponding to the (no longer 
decoupled) truncated operator H^ ^^ j^, in a spirit similar to [23]. Likewise, it is not 
essential to assume that V has common limits at xi = +00 and xi = —00; one may 
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consider also "front- type" solutions with lima;j^±oo U{xi, • ) = U±, though this in- 
troduces some additional technicalities in the analysis connected with invertibility 

4.1.4. Connection with Galerkin-based Evans functions. At this point, adopting the 
point of view of [23], we consider ICj as an operator with a matrix- valued integral 
kernel, and acting on the subspace Cj of (O; dy; i^(M; dxi)) spanned by Fourier 
modes with wave-number of modulus less than or equal to J, that is, on 



jCj = fw{x)= J2 «^m(a;i)e'™-^ WmeL''{R;dxi)\. 

^ \m\<J ^ 



(4.65) 



We let Nj denote the number of these modes. Using Lemma 4.10, one verifies that 
ICj{z) e S2(L2(M X n;dx)) and thus 

<ieh,L^(RxQ;dx){lL^{Rxn;dx) " = det2,£j(/£j - lCj{z)) (4.66) 

is well-defined. Equivalcntly, since the Fourier modes form an orthonormal basis, 
and hence the determinant is invariant under the Fourier transform, we compute 
(on the subspace L2(R;c!xi;C^-') oi L'^{R;dxi;f{Z<^-'^)) isomorphic to £j via the 
Fourier transform) instead of (4.35) the Fredholm determinant 

det2,i,2(R.d^^.cJVj)(7i2(R.rf^^.cJVj) -}Cj{z)), (4.67) 

where lCj{z) on L^(IR; dxi;C^-') is defined in terms of its integral kernel 

/_e-(Voo-Fbf-2)'/'|xi-x;i_ \ 



Fi{xi)Gjix[), x[>xi 
F2{xi)Gj{x[), xi > x[ 



T,_.^ , (4-69) 



as an operator with a single, matrix-valued semi-separable integral kernel, where 
Fk and Gk denote the Nj x Nj matrices 

mm., = {7i)^, iGk)m,, = {9l)_^, fc = l,2, (4.70) 

and GJ = {{Gk)j,m) denotes the transpose of the matrix Gk = {{Gk)m,j), k = 1,2. 

We briefly pause for a moment and recall the principal underlying idea here: The 
operator K-j acts on the space L^(lR; dxi; dy)) and leaves invariant its sub- 

space £j which, in fact, is isomorphic to (R; dxi ; Lj^^ {Q; dy)) , where (fi; dy) 
is the subspace of L'^{i};dy) spanned by the Nj exponentials {e*^'-'}|j|^j. Via 
the Fourier transform, Lj^j{Sl;dy) is isometrically isomorphic to C^-' viewed as 
a subspace in £^(Z'^^^). Indeed, if j G Z''"^ and \j\ ^ J, then the Fourier 
transform maps the element e*-' '^ of the basis of L^{fl;dy) into the sequence 
Sj = {Sm,j}meZ''-^ S £'^{Z'^~^). Indexing a basis in C^-' by means of the in- 
dices j £ Z'^^^, \j\ ^ J, we fix an isomorphism between and the subspace of 
£^(Z'^~-^) spanned by 6j, \j\ ^ J, and thus between C^^ and Lf^ ^{i}; dy). Clearly, 
via the Fourier transform, the operator fCj on (K; dxi ; Lj^^ (fi; dy)) then becomes 
the operator )Cj on L^{M.;dxi;C^''). In particular, if fCj (and hence ICj) is also 
of trace class, then the traces of ICj and /Cj are equal as will be used below in the 
proof of the second equality in (4.76). It is the operator ICj which finally is an 
operator with a semi-separable integral kernel. 
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Remark 4.12. The vector- valued case w € M" may be treated similarly, with Voo 
now a positive-definite n x n matrix, and Fk, Gk & ^nNjxnNj _ 

Our objective is to relate the truncated 2-modified Jost function (4.35) and the 
Evans function for the eigenvalue problem (4.7). We rewrite (4.7) in matrix form 

as 

(-^ + K,o/c«. +diag{bf }|^.|^^-z/c«. +W,7^vi.^0, (4.71) 

where j is the Fourier wave number, \j\ ^ J, ^ = ^{xi) is an C^' -valued function 
on M, diag { | j p } j^. | ^ ^ is a diagonal matrix of dimensions Nj x Nj , and Wj + Voo Ic"j 
is the matrix representation of some chosen truncation of the convolution operator 
V* ■ . This choice is to be followed consistently in both Galerkin-bascd and Frcdholm 
computations. Specifically, if ^ = ("0771)17711^7 for the eigenfunction -ip in (4.7), then 

Wj{xi)={Wj.m{xi))y^^^^^^j, XI GM, (4.72) 

where Wi{xi) are the Fourier coefficients oiW{xi, ■ ) = V{xi, ■ ) — Voo- Introducing 
the self-adjoint operator Hj in L^(R; dxi; C^-'), 

nj = -^ + Ko4=(M;dxi;C-.)+diag{|jf }| dom(Wj) = H^R;C''^), 

(4.73) 

we note that the asymptotic operator for the operator Hj in (4.73), viewed as a one- 
dimensional matrix- valued second-order differential operator in L^(M; dxi; C^-'), is 
given by 

72 

(4.74) 

and thus the operator 

ICjiz) = -{nf -zIr^.^^.a,^.c-j)y'Wj, z€C\a{nf), (4.75) 
is the Birman-Schwinger-type operator (cf. (4.12)) for the pair of the truncated 
operators Hj and W j ^ . 

Lemma 4.13. Let z £ C such that z ^ (t(7Yj''') . In addition, assume that Wj-m £ 
£1(L2(M; dxi)) n C(R) for all m,j £ Z, \m\, \j\ J. Then the operator Z.j{z) is of 
trace class on L^{R;dxi;C^-'), the operator K,j{z) is of trace class on L^{Rxfl;dx), 
their traces are equal and given by the following expression Qj{z): 

Qj{z) = tTL2^^^Q.ax){ICj{z)) = il^LHR;dxuC'^i){^ji^)) (4-76) 

= -U [ dxWix)) ^(K» + |j|2_^)-i/2. (4.77) 

^ \ JrxQ J I .|^ 7 

Proof. By Lemma 4.10 {ii), fCj{z) and each K,m,j{z) in (4.47) is a trace class opera- 
tor on L^(]R X f2; dx). In addition, each of the integral kernels of ICj{z) and ICm,j{z) 
is continuous on the diagonal. Thus, [2, Corollary 3.2] applies and hence 

trL2(Rxn;da:)(^j(2)) = ^ ^'^L'^{MxQ;dx){l'^m,j{z)) 
\m\,\j\^J 
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|m|,bKJ' 



proving (4.77). 

Finally, denote by Pj the orthogonal projection in L^(]R x Q;dx) onto £j and 
by Qj = lL^{RxQ;dx) — 'Pj the complementary projection. Since £j is a reducing 
subspace for ICj, 

ICjPj = VjKj, (4.79) 
one can write Kj in L^(M x Q; da;) in terms of the 2x2 block decomposition 

f^A CO ^^z) 0\ 
Kj{z)= [ ''1-"^^^)^ ' \ (4.80) 



with respect to the decomposition 

i^(K X O; rfa;) = Pji^(K x f7; rfa;) © QjL'^{^ x fi; da;) 



(4.81) 



Since ^jl^ancpj) is unitarily equivalent to Kj{z) via the Fourier transform, (4.80) 

implies that JC.j{z) and K.j{z) arc trace class operators at the same time, and it 
also implies equality of the following traces: 



trz,2(Kxa;dx)(^j(^)) = *r£j(^j|ranCPj)(^) 



(4-82) 

proving the second equality in (4.76). □ 

Remark 4.14. We emphasize that the sequence {Qj{z)}j^o diverges as J ^ oo 

for d ^ 2. The latter fact docs not permit us to use in the subsequent analysis the 
non-modified Jost function J^j{z) = deti,2(Rxn;da:)(-fL2(Rxn;dx) - l^j{z)) and pass 
in (4.93) to the limit as J — > oo; instead, it forces us to work with the 2-modified 
Jost function J^2,j{z) defined in (4.35). 

Denote by £j the Evans functions for the onc-dimcnsional approximate system 
(4.7) obtained by Galerkin approximation/ Fourier truncation at the level \j\ ^ 
J, for simplicity of discussion normalized as described in [11] to agree with the 
corresponding (one-dimensional) 2-modificd Frodholm determinant. Following the 
approach of [11], we recall that the Evans function £j{z) is defined as a 2Nj x 2Nj 
Wronskian 

£j{z) = det{y) = det{y+, . . . , y^^^), (4.83) 
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where the {2Nj x l)-vectors , . . . , y2Nj bases of solutions decaying at xi = 
+00, respectively, at xi = — oc, of the first-order system equivalent to the second- 
order differential equation (4.71), lying in appropriately prescribed directions at 
spatial infinity. The solutions are chosen in [11] in a way that £j does not depend 
on the choice of coordinate system in C^^' and docs not change under similar- 
ity transformations of the system. Specifically, let us introduce the {2Nj x 2Nj) 
matrices 

^=(^ V)' ^(-)=U(.0 2)' (4.84) 
where, for brevity, we denote 

h{z) = {Voo - z)/c«. +diag{|jf }|^.|^^. (4.85) 

Now dy/dx\ = (^ -|- i?(a;i))3^ is the first-order system equivalent to the second- 
order differential equation (4.71). We introduce on the space L^(M; dxi\€?^-') the 

first order differential operators G^j^ = (d/dxi) — A and Gj = (d/dxi) —A — R{xi) 
and the respective Birman-Schwingcr-typc integral operator 

jCj{z) = -{{d/dxi) - A)-^R{xi). (4.86) 

According to the main result in [11], the Evans function £j{z) is equal (up to the 
explicitly computed factor e^^ '(^\ that is non degenerate and analytic with respect 
to z) to the 2-modified Frcdholm determinant of the operator I (^i-dxi-c^'^' j ) ^ i/^) 
on L"^ (M; dxi ; C^'^-' ) corresponding to the first-order system mentioned above. As 
we will see next (cf. also [13, Theorem 4.7]), this 2-modified Predholm determinant is 
equal to the 2-modified Predholm determinant of the operator lL'^(R-dxi;C'^j) ~^{^) 
on I/^(M; dxi;C^-') corresponding to the second-order operator Hj in (4.71). Thus, 
we obtain evidently that the Evans function £j coincides with the following (non- 
modified(!)) Predholm determinant (that is, with the Jost function): 

J^j{z) = deti2(K;<jxi;C'Vj) (42(M;dxi;C''J) + (^J ^ " zh'^iR;dxi;C^J)) ^V^j) ■ (4-87) 

Thus, we have the following main result, extending the one- dimensional theory 
of [11]. One of its main points can be explained as follows: Zeros of the Evans 
function £j or, equivalently, of the Jost function J^j, are the eigenvalues of the 
operator Hj. The eigenvalues are also zeros of the 2-modified Predholm determinant 
J^2..j- The modified and nonmodified determinants are equal up to the nonzero 
exponential factor e^-' , where &j is the trace described in (4.77). From this point 
of view the use of the nonmodified determinant J^j (or £j) and the 2-modified 
determinant J^2,j are equivalent, as far as finding the eigenvalues of Tij is concerned. 
However, the nonmodified determinants have an advantage because the sequence 
{^2,j} converges to .F2 as J — > 00, while the sequences {^j}, {Sj}, and {9j}, 
all diverge. Thus, for the truncated problem, the use of the 2-modified Fredholm 
determinants appears to be more natural than the use of the Evans function. 

Theorem 4.15. Let z G C such that z i cr{nf^) and assume that Wj—m G 
e^{L'^{R;dxi))r]C{R) for allm,j G Z, |to|, |j| < J. Then the Galerkin-based Evans 

function £j{z), the Jost function Tj{z) = dctL2(jg^xQ.cix){lL^Kxn;dx) and 
the approximate modified Fredholm determinants for the operators in (4.75) and 
(4.86) are related as follows: 

^2,j{z) = det2,L^{lS.xn;dx){lL^(RxQ;dx) - (4.88) 
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= det2^L'^(K:dxi;C^"j) {lL^(S.;dxi;C^"j) " ^j(^)) (4.89) 

= e®^^'^£jiz) (4.90) 

= e®-'^'^'^detL2(^R^n;dx)ilLHRxn;dx) - (4.91) 

= e®^(^)^j(0), (4.92) 

where &j{z) is the trace of the operator )Cj{z) given in formula (4.77). 

Proof. Since )Cj{z) is of trace class on L^(R x f2; dx), Qj{z) is the trace of ICj{z), 
and Tj{z) is just a notation for deti2(Rxn;dx)(-?^L2(Rxn;da;) - the equality of 

(4.88), (4.91) and (4.92) trivially follows from (cf. (2.41)) 

dei2,L'^(Rxn;dx){lL'^[WLy.(l;dx) " 

(4 93) 

= deti2(Kxf,^rf,)(/i2(Rxj^^d,) - /C,7(z))e'^-^(«x":-)('= .'(-)). 

To show that the modified Fredholm determinants in (4.88) and (4.89) are equal, 
we will utilize an idea from the proof of [22, Proposition 8.1] (see also a related 
result in [13, Theorem 4.7]). We introduce the following operator matrices acting 

on L'^{^-dxi;'C'^^'): 



/ I J ' ^ 0^ 

where / = lL^(]g^-dxi-C'^j) h is defined in (4.85). In addition, we use the related 
similarity transformation in (4.84) and (4.86) to define the following matrices and 
operators 

^(1) = TAT-\ R^^\xi) = TR{xi)T-\ K}-}\z) = TKj{z)T-'^. (4.96) 
A short calculation reveals: 

=diag{-/ii/2,/iV2}, (4.97) 
R^^\xi) = (4.98) 

Changing the order of multiplication of the operators blocks in (4.99), one infers 

Thus, using the standard determinant property (2.2) and recalling (4.75), we con- 
clude that (4.88) and (4.89) are equal, 

= det2^i^2(^.dxi;C^"J){h'^(R;dxi;C^^J) ~ ^J^(^)) 

= det2,L2(^.dxi;C^J){lL^(R;dxi;C^J) " {{d^/d^xl) - h)~'^Wj) 

= det2,L2iM-dx,;C-j){l + {nf^ - ziy'Wj) 

= det2^i2(R.d2;j.cNj)(-'"l,2(R.d2;j.cNj) " }Cj{z)) 
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= det2,L2(RxO;d:r)(-^L2(RxJ2:dri;) - ICj{z)). (4.101) 

Finally, to show that (4.89) and (4.90) are equal, we apply the similarity trans- 
formation (4.96) and replace the differential equation dy/dx\ = {A + R{x\))y by 
dy/dxi = (^(i)+i?(i)(.Ti))3^. We will now use one of the main results of [11]. Since 
the real part of the spectrum of h^/"^ (for h defined in (4.85)) is positive due to our 
convention Im(.2) > 0, the unperturbed equation dy/dxi = A'^^^y, due to (4.97), 
has the exponential dichotomy on R with the dichotomy projection Q defined in 
(4.95). Therefore, according to [11, Theorem 8.37], under assumption (4.13) we 
have the formula 

deh,LHR-M^-J) {lL-(R;dx,;&-J) ' {z)) = 6^'^''^ £j{z), (4.102) 

where ^j{z) is defined as follows: 

^ poo nO 

ej{z) = / dxi trc2iv, [QR^^\xi)] - / dxi trc2«^ [(/c2«^ - Q)R^'\xi)] . 

Jo J-oo 

(4.103) 

Using (4.95) and (4.98), it follows that Qj{z) — Qj{z), completing the proof. □ 

Remark 4.16. We have a similar result in the vector- valued case mentioned in 
Remark 4.12 also in the front-type setting discussed in Remark 4.11. 

4.1.5. Alternative computation. For its own interest, and for reference in the fol- 
lowing subsections, we mention an alternative method of computing .7^2, j directly 
from the reduction of [13], where the Jost function has been computed, without 
carrying out the full analysis of [11] relating this to the Evans function. Comparing 
(4.69), and [13, (1.17)] with a = 1, we have the following representation: 

^2,j{z) = det2^c"j (^^C^J - J dxi dx[ G2(xi)(/c2Nj -|- J{xi,x[))F2{x[)^ , 

(4.104) 

where 

J{xux[) = C{xi)U{xi)-'^U{x[)B{x[), (4.105) 
B={GJ -Gjf, C={F^ F2), (4.106) 

and U is any nonsingular solution of the first-order system 



dW(x 



^ = A{xi)U{xi). (4.108) 



dxi 

The formulation (4.108) is not numerically useful, since the off-diagonal ele- 
ments of A are exponentially growing with rate of order e'^l^^^l. However, noting 
that Gj2i^A)^ respectively, Fi, 2(2^1) factor as diag {e'F(^~+l-''l^~^)^''^=^'iyVj}, respec- 
tively, diag{(yoc + jjf - z)-i/2e±(^<»+bf we may reduce (4.108) by the 
coordinate change V = DU, with 

''='^^n(Kc+b?-.)v- (Ko+iip-.)^^}' ^'-'"'^ 
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to a bounded-coefficient system (cf. (4.85)) 



dVjxi) 
dxi 



Ab{xi)V{xi), 

fh^/^ + h-^/^Wj{xi) h-^/^yvj{xi) 

\ -h-^/^Wj{xi) -/iV2 _ h-^/^Wj{xi) 



(4.110) 



Ab{xi) 



of a form readily solved by the same techniques used to solve the first-order eigen- 
value ODE for basis solutions in (4.83). Indeed, this can be recognized as 
essentially the same ODE. 

Remark 4.17. Likewise, one might start with the Evans formulation (4.71), writ- 
ten as a first order system, and try to precondition by factoring out the expected 
asymptotic behavior, to obtain essentially system (4.108). That is, the operations of 
preconditioning (viewing the Fredholm formulation as an analogous precondition- 
ing step of factoring out expected spatially-asypmptotic behavior) and reduction 
to ODE essentially commute, at least in this simple case. 

4.1.6. Stability index computation. Following the approach of Section 4.1.5, compu- 
tation of the multi-dimensional stability index can be carried out in the same way, 
with no additional complications. For. exactly as in (3.64) of the one-dimensional 
case (but using (2.52)), one has the formula 

•^2(0) = det2,i,2(][{xQ;dx)(-^L2(Mxa;dx) "^0 - Pq) detpg2:,2(][txQ;da;) (-Po^i-Pq), (4-111) 



with fCo = A^(0) = — -fffip uv and JCi = /C'(z)|^_q = — -f^Qp uv. (One can 
remove the closure symbols in the last two expressions since all operators involved 
are bounded.) 

The second inner-product-type factor is straightforward to evaluate, requiring 
only approximation of the eigenfunctions of the operator 7 — /Co and its adjoint 
corresponding to the zero eigenvalue, which in many cases are known from the 
outset. We recall that in the present case, the eigenfunction $ of the operator 
/ — /Co is dU/dxi; here U is the standing wave, see Section 4.1.1. The eigenfunction 
$ of the operator (7 — /Cq)* may be deduced by the fact that the original differential 
operator L is self-adjoint; specifically, ^7 — (Hq^^) ~ Voo)^^ = implies, by 

self-adjointness of (l + {H^l)~^) and {V - V^), that $ = 77^°^$ is indeed the 



The first factor in (4.111) on the other hand is the characteristic determinant of 
a rank-one perturbation at z = 0, so can be approximated as in Section 4.1.5 using 
Galerkin approximation/semi-separable reduction by a finite dimensional determi- 
nant. Precisely, combining the steps of Sections 4.1.2 and 4.1.5, one reduces the 
computation at the J-th Galerkin level to the evaluation of a 2{Nj + 1) x 2{Nj+ 1) 
determinant, obtained by solving a 2{Nj + 1) x 2{Nj + 1) ODE system 




dl^{xi) 
dxi 



A{xi)U{xi), U{xi) G c2(^^+i)x2(^^+i), 



(4.112) 
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where, similarly to (4.108), 

[-GJF, -GJfJ' (4.113) 

and $ is the eigenfunction of the operator JCq corresponding to the eigenvalue 1, 
that is, )Co^ = ^- Making the change of coordinates U = DV, D = diag{i), 1}, see 
(4.109), we obtain a system dV{xi)/dxi = Ai,{xi)V{xi) with bounded coefficient 
matrix Af, that can be numerically solved by standard techniques used to compute 
the Evans function. 

By comparison, if one follows the existing Galcrkin methods, working with an 
approximate truncated system at level J, one must face the difficulty that zero 
eigenvalues for the exact system perturb to small but in general nonzero eigen- 
values of the approximate system, making difficult a straightforward numerical 
computation without further analytical preparations. On the other hand, the usual 
analytic preparations (see [1],[30]) involve solving variational equations about the 
zcro-cncrgy eigenfunction and also projecting out the zero cigcnmodcs from other 
modes to obtain a well-conditioned basis. These do not appear to be real obstruc- 
tions to the computation, but are at least complications. Perhaps for this reason, to 
our knowledge no such computation has so far been carried out, or even proposed 
in fuU detail.^ 

The formulation of the above multi-dimensional stability index algorithm we 

thus view as a useful practical contribution of the present work, and its numerical 
realization as an important direction for further investigation. 

4.1.7. Numerical conditioning. Last, we examine the question of numerical condi- 
tioning. By Theorem 4.15, one way to compute the approximate Fredholm deter- 
minants J-^2.j is to carry out a standard Evans fimction computation as in [23]. 
However, from a numerical perspective, this might be missing the point. For, note 
that the principal, constant-coefficient diagonal, part of the coefficient matrix of 
(4.71) has entries leading to spatial growth rates ±|j| of order up to J. Com- 
puting the Evans function thus involves solution of an ODE that becomes infinitely 
stiff as J ^ 00. 

We suggest as a possible alternative, discretizing the kernel ICj in variables xi, 
x[ and directly evaluating the determinant of the resulting MNj x MNj matrix, 
where M is the number of mesh points in the x\ {x'l) discretization required to 
give a desired error bound. Noting that Ad is essentially dimension-independent for 
simple first-order quadrature (since matrix norm |W,/| is bounded, by Parseval's 

(y -|-|j|2_z)i/2 — I is of order one), we 
see that there should be a break-even point at which the cost of order (MNj)^ of 
evaluating the discretized determinant should be better than the cost ~ N{J)Nj 
of evaluating the Evans function, or, equivalently, 

N{J) > M^, (4.114) 

where N{J) denotes the number of mesh points required to evaluate the 2Nj x 2Nj 
Evans ODE to the same tolerance. For a first-order A-stable scheme, N{J) ~ M 



'However, see the interesting analysis [28] in the somewhat different spatially periodic case. 
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(note: this is determined by truncation error, which is estimated as proportional to 
second derivative of the solution as the square of the norm of the largest eigenvalues 
±J), yielding break-even at J ~ M, where typical values of M are of order ~ 100 
[18]. 

Though hardly conclusive, this rough calculation suggests at least that direct 
Fredholm computation is worthy of further study; we recall that J ~ 100 is the 
order studied in [23]. Alternatively, one might compute the Evans function not 
by shooting, but by continuation- type algorithms as suggested by Sandstede [34], 
viewing the eigenvalue equation as a two-point boundary-value problem, avoiding 
stiffness by another route; however, so far as we know, such a scheme has not yet 
been implemented. See [18] for further discussion of this approach. 

4.2. Functions with radial limits. Finally, consider standing-wave solutions U 
of (4.1) on the whole space M"^, possessing radial limits in the following sense: 
Introduce spherical coordinates x = {r,uj), r > 0, lo G S'^~^, and let U{R,u)), be 
the restriction of U to the sphere of radius R. Then considered as a function of the 
angle oj € U{R, ■ ), has an ^^(S"'"^; da;d_i)-limit as R ^ +oo. 

Remark 4.18. Assuming the hypotheses of Subsection 4.1, there exist radially 
symmetric solutions U{\ ■ |), where U is the solution of the corresponding one- 
dimensional problem with the second derivative replaced by the spherical Laplacian. 
Linearizing about U and expanding in spherical harmonics, one obtains a decoupled 
family of one-dimensional eigenvalue problems, similarly as in Remark 4.1, each 
of which possess a well-defined Evans function and stability index. In this case, 
there is no zero-eigenvalue at the zeroth harmonic (constant function), but there 
is a zero-eigenvalue of order d at the level of the first harmonic, with associated 
eigenfunctions dU{\ ■ \)/d,Xj, j = l,....d, corresponding to translation- invariance 
of the underlying equations. These may be treated similarly as in Remark 4.1. 
(However, we note that this involves an Evans function on the half-line [0, +oo), 
which involves some modifications and will be analyzed elsewhere.) 

Remark 4.19. In the general case, Galerkin approximation in spherical harmonics 
yields a finite-dimensional system for which an Evans function and stability can 
again be defined, similarly as in Remark 4.3. 

4.2.1. Fredholm determinant version: radial case. In the simplest situation that U 
has a single limit as — > -|-oo, the operator Hq^^ is again constant-coefScient, and 

the procedure of Section 4.1.2 leads again to expansion of H^^], ™ ^ coimtable sum 
of operators with semi-separable integral kernels corresponding to the restrictions 
to different spherical harmonics. In the general case, we may proceed instead by 
Galerkin approximation as described in Remark 4.11. 

Remark 4.20. The common feature of the problems discussed is the presence of a 
single unbounded spatial dimension (axial for cylindrical case, radial for the radial 
case), along which the semi-separable reduction is performed. In principle, one 
could treat still more general problems by truncation/disrctization of a continuous 
Fourier integral. In this setting, the reference to a concrete object in the form 
of a Fredholm determinant might become still more useful for numerical valida- 
tion/conditioning, as compared to ad hoc constructions like those in Remarks 4.1 
and 4.3. However, it is not clear that there would be a computation advantage to 
doing so. 
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4.3. General operators. We recall from [11], that it was necessary for general 
first-order operators to relate the Evans function and a 2-modified determinant 
already in the one-dimensional case, since the Birman-Schwinger kernel is for first- 
order operators only Hilbert-Schmidt (indeed, this is one of the key insights of [11]). 
Likewise, for more general operators involving a first-order component, in particu- 
lar those arising in the study of stability of viscous shock solutions of hyperbolic- 
parabolic conservation laws appearing in continuum mechanics [47] , it is necessary 
in dimensions d = 2 and 3 to relate the Evans function to a higher modified Fred- 
holm determinant, since the Birman-Schwinger kernel is no longer Hilbert-Schmidt. 
Flow in a cylindrical duct has been studied for viscous shock and detonation waves 
in [40], [41]. 
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